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LIST  OF  SYMBOLS 

t Geometrical  Quantities 

r - the  space  vector  in  the  rectangular  coordinates 

Xi,  x2  and  x3  • the  components  of  the  space  vector  in  the 
rectangular  coordinates 

x - the  coordinate  along  the  electron  flow 

t -•  the  time  variable 

z - cylindrical  coordinate  along  with  the  electron  flow 
r - the  radial  distance  from  the  z-axis 

R - the  radial  distance  in  the  spherical  coordinates  which  takes  the 
center  of  the  anode  orifice  as  the  origin 

0 - the  polar  angle  of  the  spherical  coordinates 

$ ••  the  azimuthal  angle  of  the  spherical  coordinates 

6 - an  incremental  polar  angle  from  0 

Cl  - the  solid  angle  in  the  unprimed  coordinates 

Cl'  the  solid  angle  in  the  primed  coordinates 

la  - the  length  of  the  cylinder  to  the  left  of  the  anode  plane 

12  - the  length  of  the  cylinder  to  the  right  of  the  anode  plane 

D - the  cathode-anode  spacing  of  the  diode 

01  the  polar  angle  that  the  expanded  beam  :at  the  back  wall  (due 
to  electron  lens  effect  alone)  would  subtend  at  the  origin 

02  - the  polar  angle  that  the  expanded  beam  at  the  back  wall  (due 

to  space  charge  effect  alone)  would  subtend  at  the  origin 

y - r/ro,  the  ratio  of  the  radius  of  the  electron  beam  at  the  z-plane 
and  that  of  the  anode  orifice. 

2 Electrical  Quantities 

V - the  potential  in  volts 

p - the  charge  density  in  the  electron  flow 

J - the  electric  current  density 
“ ♦ 

F the  force  per  unit  mass  of  electrons 

n the  number  density  of  the  electrons  at  a point  in  the  flow 
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E - the  applied  electric  field 

Ei  the  electric  field  at  the  left  end  of  the  cylinder 

E2  - the  electric  field  at  the  right  end  of  the  cylinder 

Er  - the  radial  components  of  the  electrical  field  across  the 
mantel  surface 

- the  anode  current  density  in  ^la/cm2 

10  the  total  current  through  the  anode  orifice 

Ip  - the  primary  electron  current  to  the  inner  wall  of  the  anode  orifice 

Is  - the  secondary  electron  current  from  the  inner  wall  of  the  anode  orifice 

is  - the  secondary  electron  current  per  unit  area  of  plane  surface 

j (8,6)  - the  current  between  the  concentric  cones  with  the  apex  angles  20 
and  2(0  + 6)  per  unit  area  of  the  anode  orifice 

N 

J(0,6)  - the  current  density  at  the  Z-plane  between  the  polar  angles  0 and 

(0  + 5) 

J ( 0)  - the  current  density  at  the  Z-plane  with  a polar  angle  0 

Js(6)  - the  current  density  at  the  Z-plane  with  a polar  angle  G due  to 
secondary  electrons  alone 

W - the  work  function  of  the  emission  surface 

T - the  cathode  temperature 

Rs  - the  ratio  of  the  space  charge  current  to  the  current  emission  of  the 
cathode 

3 Dynamic  Variables 

v0  - the  initial  velocity  of  the  electrons  at  the  cathode 

vx  - the  velocity  of  the  electrons  at  point  x 

H ••  the  velocity  vector  of  an  electron  in  the  rectangular  coordinates 

Ci,  c2  and  ca  the  velocity  components  with  respect  to  the  rectangular 
coordinates 

cr  - the  r-component  velocity  of  the  electron 
cz  - the  z- component  velocity  of  the  electron 
u - the  vector  of  the  peculiar  velocity  of  an  electron 
c the  mean  velocity  of  the  electrons  at  a point 
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the  mean  velocity  along  the  electron  flow 

the  summational  invariants 

the  mean  transit  time 

the  relaxation  time 

the  mean  free  path 

the  electron  temperature  generally  a tensor  at  a point  in 
the  electron  flow 

02.  82  - the  diagonal  components  of  the  electron  temperature 
the  z-component  of  the  electron  temperature 
the  r-component  of  the  electron  temperature 

- the  velocity  distribution  functions  of  the  electrons  at  a 
point  in  the  electron  flow 
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4.  Constants 

m - the  mass  of  an  electron  (9.10  (10)-2e  gm.) 
e - the  charge  of  an  electron  (4.80  (10)-1°  e.s.u.) 
k - a Boltzmann  constant  (1.38  (10)~le  erg/degree) 
r0  - the  radius  of  the  anode  orifice  (0.00675  inch) 

Z - the  coordinate  of  the  back  wall  of  the  anode  chamber,  (Z  = 1.200  inch) 

tty  - the  thickness  of  the  anode  wall  (r0/3) 

Q - the  proportional  constant  in  the  Chi Id- Langmuir  space-charge  law 
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r\  _ -fflU—  , the  ratio  of  the  mean  energy  and  the  internal  energy  of  the 

2k0r  electrons  at  a point  in  the  electron  flow 

e - the  normal-energy  constant  of  an  emission  surface 

A,  B,  G,  L and  M - constants  in  the  equations 

h,  h'  and  H - the  proportionality  constants  in  the  equations  of  secondary 
electron  emission. 
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PART  I 

THERMODYNAMICS  OF  THE  ELECTRON  FLOW 

1 


1 . INTRODUCTION 


In  the  following  pages,  the  electron  flow  in  a diodic  arrangement 
will  be  treated  from  a general  point  of  view.  To  elucidate  an  essential 
feature  of  the  physical  behavior  of  the  diode,  a different  viewpoint 
will  be  taken  at  the  very  beginning.  In  contrast  to  the  usual  represen- 
tation where,  at  a constant  cathode  temperature,  the  anode  potential  is 
varied  and  the  anode  current  is  taken,  we  shall  consider  here  the  case 
where  the  anode  potential  is  kept  constant  but  the  cathode  temperature 
is  raised.  Thus,  if  a constant  positive  anode  potential  is  applied  and 
the  cathode  temperature  is  gradually  raised,  the  current  in  the  diode 
will  first  ascend  quickly  and  then  reach  saturation  with  further  in- 
crease of  the  cathode  temperature  This  fact  leads  to  the  conclusion 
that  many  of  the  emitted  electrons  are  forced  to  return  to  the  cathode 
when  the  current  reaches  the  saturation.  The  interpretation  of  this 
behavior  is  that  in  spite  of  the  existence  of  a positive  potential  on 
the  anode,  the  high  accumulation  of  space  charge  in  front  of  the  cathode 
may  still  depress  the  potential  there  so  low  that  a considerable  amount 
of  emitted  electrons  will  return  to  the  cathode. 

In  the  process  of  developing  an  adequate  quantitative  description 
of  the  physical  behavior  of  the  diode,  two  phases  can  be  distinguished 

in  the  papers  of  the  early  workers  in  this  field  which  successfully 

explain  several  features  of  the  diode.  The  problem  was  first  studied1 
under  the  assumption  that  all  the  electrons  leave  the  cathode  with  the 
same  initial  velocity  and  are  accelerated  in  exactly  the  same  manner  in 
the  diode.  In  other  words,  the  electron  velocity  at  any  point  in  tne 
diode  is  considered  as  being  single- valued  Then  the  problem  may  be 
attacked  only  from  Poisson's  equation,  the  energy  equation  and  the  re- 
lation between  current,  charge  and  velocity,  namely 

_ 4np  (1) 

dx2 

1 2 1 2 _ v 

2 mvx  _ 2 mV°  “ eV  (2) 

J - pvx  (0) 

where  x is  the  distance  from  the  cathode  measured  perpendicular  to  the 
plane  of  the  cathode,  vx  is  the  velocity  of  the  electron  at  a point  x, 
vQ  is  the  initial  velocity  of  the  electron  at  the  cathode,  V is  the 
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potential  at  any  point  determined  with  respect  to  the  cathode,  e and  m 
are  the  charge  and  mass  of  the  electron,  p and  J are  the  charge  and 
current  densities.  The  solution  thus  obtained  is  certainly  not  general 
and,  at  best,  rnay  be  applicable  only  to  the  case  where  all  the  electrons 
at  any  one  point  are  moving  in  the  same  direction. 

The  essential  work  in  the  second  phase  was  done  by  Epstein2,  Fry3 
and  Langmuir4.  Instead  of  assuming  the  electron  velocity  to  be  single- 
valued in  the  diode,  they  considered  a velocity  distribution  of  the 
emitted  electrons  at  the  cathode  surface.  The  usual  representation  of 
the  potential  distribution  in  a diode  is  given  in  Fig.  1.  The  potential 
depression  in  front  of  the  cathode  is  particularly  emphasized  in  this 
figure.  If  the  potential  minimum  is  found,  one  can  determine  precisely 
the  trajectories  of  all  electrons  after  they  are  emitted  from  the 
cathode.  The  interpretation  was  that  those  electrons  with  initial 
kinetic  energy  higher  than  the  potential  minimum  would  come  through  to 
the  anode  and  those  with  legs  initial  kinetic  energy  would  return  to  the 
cathode.  It  can  be  said  to  the  merit  of  Epstein,  that  he  brought  up  two 
questions:  a:  What  distribution  function  one  must  use  to  describe  the 
current  emission  in  terms  of  the  initial  velocity  of  the  emitted  elec- 
trons; b:  Where  one  must  consider  the  problem  of  encounters  between 
electrons  and  electrons.  Unfortunately,  he  did  not  attempt  to  solve  the 
electron  flow  problem  taking  encounters  into  consideration,  but  rather 
assumed  that  the  electrons  are  unaffected  by  each  other  throughout  the 
electron  flow  in  the  diode.  On  the  other  hand,  he  was  forced  to  treat 
the  electron  gas  between  the  two  electrodes  as  being  in  on  equilibrium 
state  which  would,  of  course,  demand  that  an  energy  exchange  take  place 
between  the  particles  of  which  this  gas  is  composed. 

To  form  a description  of  the  gas  under  consideration  which  would 
come  as  close  as  possible  to  the  actual  behavior  of  the  electrons  in  the 
diodic  arrangement,  one  should  consider  all  the  effects  expected  to 
contribute  essentially  to  the  overall  behavior  of  the  electron  gas  under 
the  conditions  in  question.  Three  microscopic  phenomena  will  be  con- 
sidered as  the  decisive  factors  in  the  determination  of  the  macroscopic 
behavior  of  the  electron  gas. 

First,  or.c  has  to  consider  the  change  of  particle  density  with 
respect  to  the  space  coordinates  due  to  the  accelerated  motion  of  the 
particles  between  the  cathode  and  the  anode.  Examination  of  the  parti- 
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FIGURE  1 CHARACTERISTIC  FORM  OF  THE  POTENTIAL  DISTRIBUTION 
IN  A PARALLEL  PLANE  DIODE 


ASSOCIATED  TRANSMISSION  COEFFICIENT 


FIGURE 
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c les  at  any  fixed  location  in  the  flow,  shows  ear'  jf  them  moving  in  a 
different  direction  with  a different  velocitv.  Toe  usual  method  is  to 
decompose  the  velocity  of  each  particle  into  two  vectors,  one  being  the 
mean  velocity  which  is  the  same  for  all  the  particles  at  that  location, 
and  the  other  being  the  spread  velocity  of  that  particular  particle  with 
respect  to  the  mean  velocity.  Thus  the  particles  at  that  location  in 
the  flow  will  have  a mean  velocity  on  which  the  velocity  spread  is 
superimposed . 

Next,  one  has  to  consider  the  external  forces  acting  on  the  parti- 
cles and  the  change  of  particle  density  with  respect  to  their  three 
velocity  components.  These  three  velocity  components  together  with 
their  corresponding  space  coordinates  constitute  in  the  usual  way  the 
phase  space  (p-space)  for  that  gas.  The  external  forces  acting  on  the 
particles  are  those  due  to  sources  outside  the  location  of  the  [articles 
in  consideration. 

Finally,  in  case  of  an  electron  gas,  the  internal  forces  acting  on 
the  particles  in  question  must  be  also  considered.  These  interanl 
forces  are  produced  by  the  interactions  between  the  particles  which  are 
very  close  together.  From  a statistical  point  of  view,  it  is  believed 
that  the  effect  of  these  interactions  can  best  be  described  by  a concept 
of  encounters  between  the  different  particles.  Since  the  treatment  of 
that  problem  in  this  paper  includes  the  effect  of  the  internal  forces 
on  the  process  of  encounters,  the  precise  meaning  of  the  potential5 
will  be  one  from  which  only  the  external  forces  are  derivable. 

The  best  tool  to  describe  such  a gas  at  any  point  within  the  flow 
is  the  Boltzmann  equation: 


3£  = _ 

9t 
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V f + 
c 


(4) 


where  the  velocity  distribution  function  f is  the  dependent  variable, 

but  F,  c,  r and  t are  independent  variables.  It  is  clear  that  ^ is  the 

ot 

m * — 4 

rate  of  change  in  time  of  .the  velocity  distribution  function,  whereas  F 

is  the  external  force  per  unit  mass,  Vr  and  Vc  denote  the  del  operators 

with  respect  to  space  coordinates  and  velocity  coordinates  respectively, 
o f 

The  quantity  written  in  Eq.  (4)  is  the  rate  of  change  of  the  veloc- 

i ty  distribution  function  due  to  the  encounters  of  the  particles,  and 
thereby  is  the  resuit  of  all  particle  encounters  at  the  location  in 


-S- 


I 


V 

I? 


! 

K 


integral  involving  the 

Boltzmann’s  equation  is 
an  mtegro-differential  equation  which  indicates  that  the  change  of  the 
velocity  distribution  function  at  any  point  within  the  flow  comes  from 
three  sources  that  are  precisely  the  three  factors  listed  above.  Since 
the  first  term  of  the  right  hand  side  of  Eq.  (4)  points  out  the  non- 
uniformity  of  the  flow  of  the  gas  under  consideration,  the  methods 
adopted  in  this  treatise  will  be  similar  to  those  developed  by  workers 
dealing  with  non-uniform  gases.  Particularly,  the  nomenclature  will 
follow  the  one  given  by  Chapman  and  Cowling6. 

It  will  be  shown  that  the  electron  flow  between  the  c .thode  and  the 
anode  can  be  completely  described  with  four  macroscopic  quantities. 
These  four  quantities  are: 

1.  External  force  field  F 

2.  Electron  number  density  n 

3 . Mean  velocity  c 

4.  Electron  temperature  0 

The  temperature  of  a gas  in  general  may  be  expressed  as  a second-order 
tensor,  which  is  an  array  of  nine  ordered  components.  In  this  paper, 
however,  the  electron  temperature  consists  of  only  three  mutually 
perpendicular  component j , 0t,  02  and  03 , which  are  calculated  from  the 
velocity  distribution  of  electrons  in  the  rectilinear  flow.  In  a molecu- 
lar gas  where  these  three  components  are  always  equal,  the  temperature 
may  be  further  degenerated  into  one  component,  thus  a scalar. 

In  a later  part  of  this  report,  two  examples  will  be  given  of  how 
the  presented  theory  can  be  applied  to  certain  flow  problems.  In 
particular  two  cases  will  be  distinguished:  the  one  in  which  the  three 

temperature  components,  03 , S2  and  03 , are  equal  to  each  oth^r  but  vary 
only  with  the  location  in  the  flow,  which  will  be  properly  called  the 
"isotropic  flow";  the  other  one,  in  which  the  temperature  component  in 
the  direction  of  flow,  0j,  will  vary  with  the  location  in  the  flow;  but 
the  other  two  components  perpendicular  to  the  flow,  02  and  03  , will  be 
constant  throughout  the  flow,  hence  different  from  0 . Accordingly  this 
case  will  be  termed  the  "anisotropic  flow". 

In  this  treatise,  the  effect  of  magnetic  fields  will  be  neglected 
This  approach  will  finally  yield  exact  solutions  for  the  steady  state. 


consideration.  Hence  the  expression  is  an 

dt 

unknown  function  f (c,  r,  t)  of  the  problem.  So 
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The  solutions  will  degenerate  to  known  expressions  obtained  from  the 
single-valued  velocity  theory  if  the  effect  of  encounters  is  neglected. 
All  equations,  except  when  explicitly  noted  otherwise,  are  written  in 
electrostatic  units,  centimeters,  grams,  and  seconds. 
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II.  THE  PROBLEM  AND  THE  METHOD  OF  SOLUTION 

In  the  preceding  section,  it  was  pointed  out  that  an  adequate  tool 

to  describe  the  behavior  of  the  electron  gas  in  the  condition  of  diodic 

flow  would  be  Boltzmann’s  integro-di  fferential  Eq.  (4).  This  would 

account  for  mutual  interaction  of  the  particles  by  the  way  of  encounters. 

From  a practical  point  of  view,  the  condition  of  the  flow  after  it  has 

attained  its  steady  state  should  deserve  first  attention  rather  than  the 

very  general  case  of  studying  the  transients  of  this  flow,  This  enables 

» 

one  to  make  a considerable  simplification  of  the  very  complex  task  of 
solving  generally  Eq.  (4),  since  it  becomes  unnecessary  to  account  for  a 
change  in  time  of  the  VDF*  at  a particular  point  in  consideration  That 
precisely  means  that 

§*  = 0 (5) 

3 1 

if  only  the  steady  state  is  considered.  Since  from  now  on  only  the 
steady-state  case  is  considered,  the  problem  should  be  to  solve  the 
steady-state  Boltzmann  equation 

0 £ 

c • vrf  ♦ F ' Vcf  - = o (6) 

for  the  VDF  at  any  point  within  the  flow.  That  precisely  means,  to  find 
uniquely  an 

f(c,  r,  f0(r0 ) ) 

which  would  satisfy  Eq . (6),  where,  because  of  the  character  of  this 
particular  equation  --  being  linear  and  of  the  first  order  --  one  has  a 
free  choice  of  the  distribution  function  at  a particular  point  in  the 
flow.  This  boundary  condition  was  formerly  introduced  as  f0(r0). 

It  is  well  known  to  workers  in  this  field  that  an  attempt  to  solve 
generally  the  Boltzmann  equation  leads  to  insurmountable  mathematical 
difficulties,  even  in  the  simplified  form  of  Eq.  (6),  where  only  the 
steady-state  case  is  considered  In  spite  of  this  fact,  several  methods 
were  devised  to  solve  this  equation  for  some  particular  cases. 

One  of  the  early  simplifying  restrictions  imposed  upon  the  problem 
was  the  condition  of  uniformity  and  steadiness  of  the  gas  under  consid- 
eration. Thisjneans,  a.)  no  external  force,  b.)  zero  mean  velocity. 
*VDF,  hereafter  stands  for  velocity  distribution  function. 
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The  solution  for  f in  this  case  is,  of  course,  the  veil  known  Maxwellian 
VDF. 

Recently,  Kihara7  studied  the  case  where  the  gas  under  considera- 
tion was  composed  of  charged  and  uncharged  particles.  In  the  case  where 
the  number  of  uncharged  particles  is  much  larger  than  the  number  of 
charged  particles,  the  collisions  among  charged  particles  can  be  con- 
sidered negligible  in  comparison  with  the  collisions  of  charged  with 
uncharged  particles.  Solutions  are  obtained  with  a further  restriction, 
namely,  that  the  drift  velocity  of  the  charged  particles  is  small  in 
comparison  with  their  thermal  velocity.  Obviously,  the  restriction 
pointed  out  above  cannot  be  imposed  upon  the  case  considered  in  this 
paper,  that  of  the  flow  of  an  electron  gas  in  a high  vacuum  diode. 
Unfortunately,  there  is  at  present  no  precise  mathematical  expression  to 

describe  adequately  the  phenomenon  of  encounters  between  electrons. 

9 f 

Even  with  an  approximate  expression  for  the  term,  , symbolizing  the 

ot 

encounters  in  the  Boltzmann  Eq.  (6),  one  still  meets  insurmountable 
mathematical  difficulties  to  find  a solution  for  f. 

No  attempt  will  be  made  in  this  paper  to  solve  for  the  VDF  under 
particular  restrictions  imposed  upon  the  problem.  The  method  suggested 
here  will  be  to  define  a class  of  VDF’ s by  stating  properly  the  physical 
properties  of  the  flow  and  then  selecting  a representative  function  of 
this  class  of  VDF' s which  may  contain  the  "true  function11  as  one  of  its 
members.  It  will  be  seen  that  the  so  selected  VDF  will  approach  any 
function  of  that  class  within  the  third  moment. 

Tn  the  course  of  describing  the  physical  properties  of  the  flow  and 
defining  in  that  way  the  class  of  VDF’ s , any  member  of  which  will  serve 
as  an  approximation  to  the  "true  function"  up  to  the  third  moment,  one 
links  the  microscopical  quantities  defining  the  VDF  with  macroscopic 
quantities  accessible  to  measurements,  as  was  mentioned  in  the  preceding 
section . 

Let  c be  the  linear  velocity  in  coordinates  fixed  with  respect  to 
an  observer  in  the  laboratory.  Then  the  electron  number  density  in  the 
neighborhood  of  a point  r in  the  flow  is  given  by 


n(r)  = / f (c,r)  dc 
where,  here  and  in  the  future,  a notation  / ■ 


indicates  an  integra- 


f 


I 


tion  carried  out  over  all  components  of  the  vector  quantity  a in  the 
entire  3 space • 

In  accordance  with  the  usual  definition  of  the  mean  value  of  an 
arbitrary  function  with  respect  to  a distribution  function  f,  the  mean 
value  c of  the  function  c,  the  velocity  of  the  electrons  in  the  flow, 
will  be  defined  as 

c = ^ J c f (c,r)  dc  (8) 

Eqs  (7)  and  (8)  will  impose  some  restrictions  as  to  the  free  choice  of 
an  arbitrary  VDF,  since  any  suggested  VDF  must  be  consistent  with  the 
conditions  defined  in  Eqs.  (7)  and  (8).  Since  these  equations  are  the 
definitions  of  the  zero  and  the  first  moments  of  the  velocity  distribu- 
tion function  f,  at  that  state  of  the  development,  it  can  be  said  that 
such  a function  f satisfying  (6),  (7)  and  (8)  would  approximate  the 
"true  function"  up  to  the  first  moment. 

In  a later  development  it  will  be  necessary  to  look  at  the  velocity 
distribution  of  che  electrons  in  the  flow  from  the  viewpoint  of  an 
observer  moving  with  the  mean  velocity  c along  with  the  electron  flow. 
Let  u be  the  linear  velocity  in  coordinates  fixed  with  respect  to  this 
observer.  Then  the  electron  number  density  in  the  neighborhood  of  a 
point  r in  the  flow  is  given  according  to  (7) 

n (r)  - / <p  ( u,  r)  du  (9) 

where  cp(u,r)  denotes  the  VDF  as  observed  by  the  moving  observer,  and  u, 
the  peculiar  velocity  of  the  electrons  in  this  frau.e  of  reference,  is 
defined  by 

u = c - c ( 10) 

Having  made  clear  the  definition  of  the  zero  and  first  moments,  we 
now  proceed  to  consider  the  second  moment  of  the  VDF.  By  the  second 
moment  of  the  VDF  is  meant  the  following  quantity: 

/ u iT  <p  (u,  r)  du 

where  the  symbol  uu  has  been  used  to  stress  the  fact  that  any  two 
components  of  the  peculiar  velocity  u can  be  combined  to  form  a second 
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order  term  with  respect  to  their  velocity.  Hence  uu  can  be  considered 
as  a tensor  with  the  nine  ordered  components: 


h 

Ul 

»1 

u2 

Ul 

u3) 

u2 

Ui 

u2 

u2 

U2 

U3 

\ 

yu3 

Ul 

U3 

u2 

U3 

/ 

u 3/ 

This  makes  it  obvious  that  the  second  moment  has  the  properties 
tensor,  the  components  of  which,  multiplied  by  a constant  factor 
conveniently  define  a new  tensor 


of  a 
_EL 
kn  ’ 


0 (r)  - / u u <p  (u,r)  du 

kn 


(11) 


which  is  generally  known  as  the  temperature  tensor  in  the  neighborhood 
of  the  point  r in  the  flow. 

Since  our  aim  is  to  obtain  results  from  the  presented  theory  which 
are  applicable  to  practical  situations,  structures  with  rectilinear 
electron  flow  (as  concentric,  coaxial,  parallel  plane  structures)  wiii 
certainly  receive  first  attention.  Restricting  oneself  to  the  general 
class  of  structures  with  rectilinear  flow,  a convenient  choice  of  the 
coordinates  can  be  made  by  placing  one  axis  parallel  and  the  other  two 
axes  perpendicular  to  the  flow.  The  axis  parallel  to  the  flow  in  the 
following  will  be  labeled  with  index  1,  the  other  two,  2 and  3.  The 
necessary  and  sufficient  condition  to  characterize  rectilinear  flow  is  to 
say  that  the  net  transport  of  matter  is  carried  out  along  one  and  only 
one  fixed  coordinate  Explicitly,  this  means  that  the  mean  velocity,  c, 
of  the  flow  is  parallel  to  the  external  force  field  F in  the  flow.  This 
immediately  implies  a simplification  of  Eq.  (8),  which  defines  the  mean 
velocity  in  all  different  directions  cf  the  flow.  With  the  particular 
choice  of  the  coordinate  system  as  described  above,  it  follows  that 


/ u,ip(u,r)  dui  - 0 , i - 2,3  (^2) 

The  choice  of  this  particular  geometry,  with  its  symmetry  proper- 
ties expressed  through  Eq . (12),  has  a most  important  bearing  on  the 
structure  of  the  temperature  tensor  0 defined  by  Eq . (11).  A glimpse  of 
the  tensor  array  for  the  u'u1  terms  indicates  immediately  that  all  the 
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non-diagonal  terms  will  vanish  and  the  temperature  tensor  0 reduces  to 


/e  i o 

e = J o e2 

yo  o 


° \ 
o 

03  / 


(13) 


Since  each  of  the  diagonal  terms  represents  the  mean  square  of  a veloc- 
ity component,  the  electron  temperature  is  merely  a measure  of  the  mean 
inetic  energy  with  respect  to  an  observer  moving  with  the  mean  velocity 
along  with  the  flow.  This  definition  of  temperature  will  be,  of 
course,  in  agreement  with  that  in  thermodynamics  when  the  particles 
under  consideration  reach  the  state  of  equilibrium. 

Having  now  properly  interpreted  the  physical  signi f icance  of 
the  second  moment,  we  are  ready  to  turn  our  attention  to  the  third 
moment  of  the  VD£.  In  accordance  with  the  usual  definition,  one  can 

'v,/" 

write  for  M3 , the  third  moment: 


M3  = J u3  cp  (u,  r)  du 

M3  can  be  interpreted  as  the  rate  of  net  flow  of  the  kinetic  energy  with 
respect  to  an  observer  moving  with  the  mean  velocity  c along  with  the 
electron  flow.  It  will  be  seen  later  in  the  development,  but  also  may 
be  worth  noting  here  that  in  practical  situations  the  average  velocity 
with  respect  to  such  a moving  observer  is  much  smaller  than  the  mean 
velocity  c of  the  electron  flow.  The  net  flow  of  the  kinetic  energy 
with  respect  to  such  a moving  observer  is  consequently  much  smaller  than 
the  flow  of  the  kinetic  energy  with  the  mean  velocity  c.  That  is 

k0  <:  < m c 2 

2 2 


Fox  this  reason,  the  third  moment  M3  can  be  considered  to  be  zero,  thus: 

J u3  cp  (u,r)  du  = 0 (14) 

It  may  be  noted  that  every  even  VDF  would  automatically  satisfy  the 
condition  expressed  in  Eq . (14).  However,  Eq.  (14)  does  not  necessarily 
restrict  VDF’ s to  be  even  functions  of  u.  This  minor  restriction 
imposed  upon  the  problem  does  not  prevent  us  from  including  a vast 
number  of  functions  in  the  class  of  possible  VDF's  and  just  enables  us 
to  procure  solutions  for  the  present  problem. 


Since  all  the  statements  so  far  regarding  the  physical  properties 
of  t.he  electron  flow  have  been  so  general  that  they  can  be  considered 
valid  whether  encounters  or' no  encounters  between  the  particles  take 
place,  we  must  now  turn  our  attention  to  the  more  specific  knowledge 
one  can  obtain  if  one  considers  that  even  if 'encounters  have  a major 
role  in  the  physics  of  the  flow,  one  is  able  to  select  certain  quanti- 
ties which  will  not  change  before  and  after  the  encounters.  Although  no 
detailed  picture  of  the  process  of  encounters  may  ever  be  obtained,  it 
is  at  once  evident  that,  for  example,  the  number  of  particles  engaged  in 
mutual  encounters  will  remain  constant  before  and  after  their  inter- 
actions. Similarly,  it  is  clear  that  the  total  momentum  and  energy  of 
the  participating  particles  will  remain  invariant  during  this  process. 
Those  quantities  remaining  unchanged  during  the  complex  process  of 
encounters  are  generally  known  as  summational  invariants.  Their  invari- 
ability gives  some  additional  bits  of  information  to  a sufficient 
description  of  the  physical  properties  of  the  flow. 

To  state  the  three  previously  indicated  conditions  in  mathematical 

9 f 

terms  one  has  to  remember  that  the  symbol  in  the  Boltzmann  Eq . (6) 

d t 

denotes  the  change  of  the  VDF  due  to  the  particle  encounters.  Thus,  the 
statement  that  the  number  of  particles  engaged  in  mutual  encounters  will 
remain  constant  before  and  after  their  interactions,  can  be  written  as 

3 £ 

/ -f-  d?  = 0 (15) 

o t 

Similarly , the  statement  that  the  total  iriOuientum  and  energ;;  of  the 
participating  particles  will  remain  unchanged  through  the  encounter 
process,  can  be  written 

3 £ 

I md  - dc  = 0 (16) 

9t 

f ir-  d5  = 0 a?) 

Z ot 

I .... 

In  general,  for  any  function  which  should  remain  unchanged  during 
encounters,  one  can  write 

/ </■'  ^ dc  = 0 (18) 

ot 

Thus,  the  summational  invariants,  d'i,  0?  and  . in  the  three  preceding 

j 

1 


L 
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cases  under  consideration,  will  assume  the  values, 


</>!  = 1 (19) 

= me  ( 20 ) 

and 

•>3  = ^ (21) 

Having  obtained  the  summational  invariants  as  well  as  the  relations 

between  the  VDF  and  macroscopic  quantities  of  the  flow,  one  can  now 
solve  the  problem  in  a straight  forward  manner.  With  the  aid  of  </e  ’ s, 
it  may  be  noted  that  the  term, 

lei 

dt 

describing  the  effect  of  particle  encounters  can  be  reduced  to  zero  if 
it  is  multiplied  by  a summational  invariant  and  integrated  with  respect 
to  c in  the  velocity  space.  Thus  multiplying  each  term  in  the  Boltzmann 
Eq  (6)  with  the  </e ’ s and  integrating  with  respect  to  c in  the  velocity 
space,  we  have 

/ '/'j  c • V^f  dc  + J ^ F ' V-f  dc  = / </o  dc 

(22) 

According  to  Eq . (18),  the  term  on  the  right-hand  side  vanishes  in  all 
cases.  By  this  method,  the  exact  expression  of 

dt 

is  therefore  not  needed,  but  the  effect  of  particle  encounters  has 
properly  been  taken  into  account. 

Inserting,  successively,  in  Eq . (22)  the  three  values  of  the 

summational  invariants,  ipi , '/'a  and  , as  given  in  Eqs.  (19),  (20)  and 
(21),  a set  of  three  integral  equations  of  the  following  form  is  ob- 
tained. 

Case  ipx  '■  /c  " Vr  f dc  + j F V,  f dc  = 0 (23) 

Case  0 1 : m/c  c • ^f  dc  + m / c F • Vcf  dc  = 0 (24) 

Case  V'  i : “ / c!  c d'e  “ / c’  ? • Vc  f de  « 0 (25) 
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With  the  aid  of  the  four  definitions  of  the  zero,  first,  second  and 
third  moments  of  the  VDF  as  given  in  Eqs.  (7),  (8),  (11)  and  (14;,  and 
after  some  mathematical  transformations  (see  Appendix  I),  the  six 
integrals  occurring  in  the  set  of  equations  above  can  be  carried  out: 

^ (nc)  = 0 (26) 

V • (kn0 ) - mn  [F  - c • (V£)]  = 0 (27) 

Lf  n©  • — * _ _ 

V ■ [(2 — ^-1)c]  + c * [V  -(bn©)]  + kn©  * (Vc)  - mnc  • [F  - c • (^c)j  = 0 

(28) 

It  may  oe  noted  in  the  set  of  the  foregoing  equations  that  "carrying 
out"  the  integrals  was  essentially  replacing  the  microscopical  quanti- 
ties, r,  c,  f(r,c)  by  a set  of  macr oscopi ca 1 1 y definable  mean  values, 
viz.,  F,  n,  c and  ©.  Their  connection  with  the  microscopical  quantities 
was  established  in  the  definition  of  the  different  moments.  It  is  in 
this  step  of  the  development  that  one  loses  the  microscopic  view  and 
gains  the  macroscopic  observation. 

The  physical  content  of  Eqs.  (26),  (27)  and  (28)  can  easily  be 
realized.  Equation  (26)  is  the  familiar  equation  of  conservation  of 
charge  and  matter.  Equation  (27)  is  known  as  the  equation  of  motion  in 
the  hydrodynamics  of  perfect  fluids.  The  third  equation,  Eq . (28), 

essentially  expressing  the  conservation  of  energy  in  the  flow,  is  the 
new  bit  of  information  which  is  obtained  by  eidopting  the  more  general 
viewpoint  of  an  existing  temperature  tensor  © in  the  flow.  The  clas- 
sical viewpoint  is  obtained  at  once,  if  the  terms  involving  0 are 
neglected.  If  © is  put  to  zero,  Eq.  (28)  reduces  to  Fq.  (27)  which  in 
losing  its  term 

V • (kr.0) 

would  precisely  become  the  classical  Eq . (2)  mentioned  in  the  intro- 

duction . 

In  other  words,  in  the  classical  case  the  problem  was  attacked  by 
solving  Poisson’s  equation  with  the  aid  of  the  energy  equation  and  the 
relation  between  current,  charge  and  voltage.  Turning  our  attention 
back  to  Eqs.  (26),  (27)  and  (28)  which  express  the  more  general  view- 
point adopted  in  this  paper,  one  may  note,  that  up  to  this  point  only 
three  equations  are  at  our  disposal,  whereas  four  macroscopic  quantities, 
F,  n,  c and  0 are  expeced  to  be  computed.  Since  Poisson’s  equation 
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correlates  two  of  the  four  macroscopic  quantities  under  consideration, 


V • F = 4rt  n,  (09) 

viz.,  F and  n,  and  is  independent  of  the  thiee  equations  mentioned,  above, 
it  can  be  used  as  the  necessary  fourth  equation  to  solve  for  the  four 
macroscopic  quantities,  F,  n c and  0 

A brief  summary  of  this  chapter  on  the  problem  and  its  solution  may 
be  appropriate.  The  problem  consists  of  solving  Boltzmann’s  equation 
for  the  steady  state  case  and  for  rectilinear  flow,  taking  an  integral 
effect  of  encounters  into  consideration.  Solving  Boltzmann’s  equation 
means,  in  the  usual  sense,  solving  for  VDF  under  the  particular  condi- 
tions of  the  flow  under  consideration.  But  with  the  aid  of  four  de- 
finitorial  equations  for  the  zero,  first,  second  and  third  moments,  a 
class  of  VDF’ s could  be  defined  which  may  contain  the  “true  function" 
and  where  any  representative  member  of  this  class  would  approach  any 
other  member  of  this  class  of  functions  up  to  the  third  moment.  With 
the  knowledge  of  the  overall  behavior  of  a group  of  particles  engaged  in 
mutual  encounters,  expressed  by  the  three  summational  invariants,  three 
partial  differential  equations  could  be  derived,  connecting  four  macro- 
scopic quantities,  F,  n,  c and  6.  With  the  additional  information  with 
respect  to  distributed  space  charges  expressed  by  Poisson’s  equation,  a 
set  of  four  mutually  independent  equations  is  obtained,  which  represents 
the  ultimate  amount  of  information  obtainable  in  describing  the  physical 
properties  of  the  electron  flow  from  this  general  viewpoint.  It  can  be 
seen  easily  that  these  four  equations  are  not  yet  sufficient  to  solve 
for  the  four  macroscopic  quantities  under  consideration.  The  reason  is, 
of  course,  that  the  four  equations  are  written  in  vector  form  and  repre- 
sent in  the  case  of  rectilinear  flow  a set  of  four  scalar  equations, 
since  F and  c in  Eq . (27)  are  parallel  to  each  other.  The  number  of 
unknowns,  however,  is  six,  since  0 has  three  components.  Thus,  two  more 
conditions  must  be  established  to  furnish  the  set  sufficient  for  de- 
termining completely  the  desired  macroscopic  quantities.  These  two 
additional  conditions  can  be  considered  as  particular  specifications  of 
the  problem.  Two  such  problems,  as  it  was  indicated  in  the  introduction, 
namely,  the  isotropic  flow  and  the  anisotropic  flow,  will  be  separately 
treated  in  the  next  section. 

Before 
solution,  i 
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concluding  this  section  on  the  problem  and  the  method  of  its 
t may  be  pertinent  at  this  time  to  give  a representative 


P 
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function  for  the  class  of  VDF’s  which  may  contain  the  "true  function"  as 
one  of  its  members.  Considering  that  any  VDF  selected  should  be  con- 
sistent with  the  Eqs.  (8),  (11)  and  (14),  which  define  the  first,  second 
and  third  moments,  or.e  may  suggest,  for  example,  a function  of  the  form, 


>(u)  - n(-^-)3/2(e1e203  )-1/!exp[-  jl&l.  + + !£_)] 

2rck  2k  ©i  02  0n 


where  n and  0’s  may  be  functions  of  r.  This  function  is  generally  known 
as  Schwarzschi Id ’ s law8  of  ellipsoidal  distribution  of  velocities. 


Ill  TWO  ESSENTIAL  TYPES  OF  FLOW 


In  the  preceding  section,  it  was  pointed  out  that  two  more  condi- 
tions are  necessary  for  the  complete  determination  of  the  cix  macro- 
scopic quantities,  F,  n,  c,  ©i,  ©2  and  ©3.  In  searching  for  the  two 
conditions  in  question,  one  may  explore  certain  specifications  of  the 
general  problem,  which  should  be  appropriate  for  the  actual  situation 
under  consideration.  For  example,  imagine  a situation  in  which  the 
group  of  particles  untler  consideration  during  their  life  time  in  the 
diode  may  have  sufficient  time  to  complete  their  mutual  interactions. 
If  it  is  possible  to  give  a reasonable  estimate  of  the  time  required 
to  complete  the  mutual  interaction,  which  may  be  called  relaxation  time 
tp,  then  this  situation  is  indeed  realized  if  the  transit  time  t'p  for 
those  particles  is  much  longer  than  the  "relaxation  time"  tp.  The 
statement  above,  that  the  group  of  electrons  is  allowed  sufficient  time 
to  "complete  their  mutual  interactions"  means,  of  course,  that  this 
group  has  reached  its  equilibrium  condition.  From  this,  one  can  immedi- 
ately infer  that  the  three  temperature  components,  0lt  02,  03 , for  this 
particular  group  of  electrons  should  be  equal  to  each  other.  Hence,  a 
situation  as  described  above  may  be  properly  defined  by: 

01  = ©2  = ©3  (31) 

Since  all  three  temperature  components  in  this  case  degenerate  to  one 
temperature  scalar,  this  type  of  flow  may  be  called  the  "isotropic  flow". 
In  our  search  for  two  more  conditions  essential  in  solving  for  the  six 
macroscopic  quantities,  F,  n,  c,  ©i,  ©2,  ©3,  one  may  note,  with  the  aid 
of  the  set  of  four  fundamental  differential  equations  (26)  - (29),  that 
the  specification  as  described  above  furnishes  exactly  the  two  necessary 
conditions  over  Eq . (31).  Thus  a necessary  and  sufficient  set  of 
equations  has  been  established  in  order  to  determine  completely  the 
macroscopic  quantities  under  consideration. 

Similarly,  imagine  a situation  in  which  the  group  of  electrons 
emitted  fiom  the  cathode  are  so  swiftly  drawn  to  the  anode  that  no 
appreciable  amount  of  particle  encounters  can  take  place  in  the  space  of 
the  diode.  This  situation  is  realized  if  the  transit  time  tj  for  those 
electrons  is  much  shorter  than  the  relaxation  time  tp  for  this  particu- 
lar group.  Obviously,  in  this  situation,  the  group  of  particles  under 
consideration  will  not  reach  their  equilibrium  condition.  However,  an 
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observer,  moving  with  the  mean  velocity  along  the  electron  flow,  should 
in  this  case  not  be  able  to  observe  any  essential  change  of  the  electron 
velocity  distribution  in  the  velocity  components  perpendicular  to  the 
axis  of  the  flow.  Thus  the  temperature  components  perpendicular  to  the 
direction  of  the  flow  become  constant  in  the  space  interval  of  consider- 
ation. Hence,  a situation  as  described  above  may  be  defined  by: 


902  _ 903 
9xi  9xi 


(32) 


Since  all  three  temperature  components  in  this  case  do  not  necessarily 
degenerate  to  one  temperature  scalar,  this  type  of  flow  may  be  called 
the  "anisotropic  flow".  One  may  note  that  the  specification  as  de- 
scribed above  again  furnishes  exactly  two  more  necessary  conditions  in 
solving  for  our  macroscopic  quantities.  Thus,  in  this  case  too,  a 
necessary  and  sufficient  set  of  equations  has  been  established  in  order 
to  determine  completely  the  macroscopic  quantities  under  consideration. 
The  two  situations  briefly  described  above,  namely  the  "isotropic  flow" 
and  the  "anisotropic  flow",  will  be  separately  discussed  in  two  later 
paragraphs . 

Now,  we  will  turn  our  attention  to  a more  precise  determination  of 
the  two  terms  already  discussed,  viz,,  the  "relaxation  time ", , t p , and 
the  "transit  time",  tp.  In  making  a reasonable  estimate  of  the  time 
required  for  a group  of  particles  to  complete  their  mutual  interactions, 
we  shall  follow  essentially  Chandrasekhar’s  extensive  wor  k9  on  the 
problem  of  encounters  in  stellar  systems.  The  analogy  with  the  electron 
gas  is  rather  obvious  if  one  considers  that  in  both  cases  the  inter- 
action forces  drop  with  the  inverse  square  of  the  mutual  distances  of 
the  particles  in  question.  Chandrasekhar  defines  a relaxation  time  for 
a star  cloud  at  that  time:  most  of  the  particles  need  to  alter,  through 
the  process  of  their  mutual  interactions,  their  original  kinetic  energy 
by  about  an  equal  amount.  Accepting  this  definition  and  converting  the 
quantities  designed  to  fit  the  stellar  case  into  quantities  appropriate 
for  the  electron  gas,  one  finds  for  the  so  determined  relaxation  time, 
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With  expression  (33)  for  the  "relaxation  time",  one  can  evaluate  another 
quantity  of  interest,  namely,  the  "mean  free  path".  It  may  seem  im- 
proper to  define  a mean  free  path  for  the  election  gas.  Nevertheless, 
it  is  instructive  to  define  a length 


X = tR  (&&)' 

11  Ttm 


which  plays  the  same  role  for  non-uniform  electron  gases  as  the  mean 
free  path  does  in  the  classical  kinetic  theory  of  gases,  where  the 
molecules  are  idealized  as  rigid  elastic  spheres.  Of  course,  this 
"mean  free  path"  is  defined  in  the  moving  coordinate  system  and  is 
applicable  whether  the  gas  as  a whole  is  actually  mo/ing  or  not. 

Having  obtained  an  expression  of  the  "relaxation  time"  for  the 
group  of  particles  under  consideration,  one  may  now  determine  the  trans- 
it time  for  this  group  of  particle’s.  Since  the  transit  time  for  a group 
of  particles  in  a diode  is  nothing  more  than  tie  mean  life  time  of  this 
group  in  a spatial  interval  an  appropriate  definition  of  what  we  may 
call  a "mean  transit  time"  may  be  given  by 


/-b  dxi 

t t = J 

0 c (Xl) 


for  a spatial  interval  between  points  a and  b in  which  the  electron  gas 
is  under  consideration.  Since  the  mean  velocity,  c,  is  a single-valued 
function  of  Xi,  formula  (35)  is  applicable  to  the  entire  region  of  the 
aiodic  flow.  Without  the  definition  of  the  mean  velocity  (Eq.  (8)),  one 
cannot  define  the  transit  time  of  electrons  in  all  cases,  especially  of 
the  emitted  electrons  which  are  eventually  returned  to  the  cathode. 

The  fundamental  distinction  between  the  two  situations  described 
above  --  the  isotropic  and  the  anisotropic  flow  --  is  that  for  isotrop- 
ism,  the  relaxation  time  (tR)is  is  very  short  in  comparison  with  the 


transit  time  t-i 


(tph 


whereas  for  anisotropism,  the  relaxation  time  (to)  should  be  much 

Q n 


longer  than  the  transit  time  ty 


This  allows  one  to  draw  an  important  conclusion  with  respect  to  the 
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number  density  in  both  cases,  since  tp  is  about  inversely  proportional 
to  the  number  density  of  the  electrons  under  consideration  (see  Eq. 
(33)).  Assuming  the  transit  time  ty  for  both  cases  to  be  roughly  the 
same,  the  two  cases  may  also  be  distinguished  by  the  inequality 

(to)  > > > (tp) . 

or  using  Eq.  (33),  which  expresses  the  relaxation  time  in  terms  of  the 
parameters  describing  the  electron  gas  under  consideration,  one  finds 

(n)an  * < < <n>is 

This  strong  inequality  suggests  that  Poisson’s  equation  will  degenerate 
to  Laplace’s  equation  in  the  case  of  the  anisotropic  flow. 

Tile  specifications  given  for  the  two  situations  are  based  upon  the 
determination  of  the  "relaxation  time"  and  the  "mean  transit  time"  of 
the  group  of  particles  in  question.  In  concluding  the  introductory  part 
of  this  section  on  the  two  essential  types  of  flow,  it  might  be  well  to 
list  all  the  equations  which  will  be  used  in  solving  the  problem- 
However,  for  the  sake  of  simplicity,  we  shall  confine  ourselves  to  the 
parallel  plane  diode.  For  the  electron  flow  in  the  parallel  plane  diode 
where  the  variations  of  the  macroscopic  quantities,  F,  n,  c,  0i’s  are 
only  in  the  direction  of  the  flow,  the  general  equations  derived  for  the 
electron  flow  are  reduced  to  simpler  forms.  The  general  equations 
(26) -(29)  expressed  in  their  reduced  form,  together  with  the  conditions 
of  the  specifications  as  mentioned  above,  are  listed  as  follows: 


(nv)  = 0 (36) 

dx 

d-  (kn0a)  - mn  (F  - v dx)  = 0 (37) 

dx  dx 

J_[k£y.(3  @1  + 02  + 03)]  - TOY  (F  - V dx)  =0  (38) 

dx  2 dx 

In  the  case  of  the  "isotropic  flow", 

F = 0 (39a) 

dx 

-d-  Gj  = -d-  03  = 0 (40a) 

dx  dx 
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In  the  case  of  the  "isotropic  flow", 


flow", 

-d-  F = 4tc  n 

dx  m 

(39b) 

01  = 02  = 03 

(40b) 

where  F,  v and  x are  used  instead  of  F,  c and  xt. 

The  problem  now  consists  of  solving  simultaneously  the  set  of  the 
above  equations,  (36)  - (40).  For  the  electron  flow,  Eqs.  (36)  - (38) 
are  valid  in  general.  In  the  case  of  the  anisotropic  flow,  Eqs.  (39a) 
and  (40a)  are  added,  whereas  in  the  case  of  the  isotropic  flow,  Eqs. 
(39b)  and  (40b)  are  added.  However,  some  insight  can  be  gained  immedi- 
ately by  integrating  Eq  (36)  with  respect  to  x.  Equation  (36)  now 
becomes 

n v = j ( 41 ) 


where  j is  an  integration  constant  of  the  equation.  Since  the  displace- 
ment current  is  zero  in  the  steady-state  case,  je  becomes  the  actual 
electric  current  density  in  the  diode  circuit.  Equation  (41)  gives  the 
functional  relationship  between  the  number  density  n and  the  mean 
velocity  v of  the  electrons  in  the  flow. 

Another  possibility  of  establishing  a functional  relationship 
valuable  foi  both  cases  would  be  to  eliminate  the  external  force  F from 
Eqs.  (37)  and  (38).  Eliminating  F from  these  equations  and  making  use 
of  Eq.  (36),  one  obtains 


0 1 dx 


(0i  +■  0?  + 0S) 


+ 2.  dx  = q 

v dx 


(42) 


In  the  case  of  the  "anisotropic  flow"  where  the  temperature  component^, 
02  and  03,  are  independent  of  x,  Eq.  (42)  yields 

0i  (x)  [v(x)]2  = 0i  (o)  fv(o)]  2 = constant  (43) 

where  0i(o)  and  v(o)  are  the  x-component  of  the  electron  temperature 
and  the  mean  velocity  at  the  emission  surface.  The  macroscopic  quanti- 
ties at  the  emission  surface  are  known  as  boundary  values,  the  evalua- 
tion of  which  will  be  given  later.  Similarly,  in  the  case  of  the 
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"isotropic  flow"  where  the  components  of  the  electron  temperature  are 
equal  to  each  other,  Eq.  (42)  then  yields 

[®i(x)]3  [v(x)]  2 = [6i(o)]J  [v(o)]  2 = constant  (44) 

These  two  equations  indica  u S vliH  t the  electron  temperature  in  the 
electron  flow  of  a parallel  plane  diode  decreases  as  the  mean  velocity 
of  the  electrons  increases  (Appendix  II).  If  one  could  find  the  mean 
velocity  as  a function  of  x in  the  electron  flow  under  consideration, 
the  variation  of  the  electron  temperature  along  with  the  electron  flow 
could  be  determined  from  Eqs . (43)  and  (44). 

In  the  following,  we  shall  consider  separately  the  two  types  of 
flow;  first,  the  "anisotropic  flow"  and  second,  the  "isotropic  flow". 

3. 1 Anisotropic  Flow 

As  mentioned  before,  the  problem  now  is  to  solve  simultaneously 
Eqs.  (36),  (37),  (38),  (39a),  and  (40a)  for  the  case  of  the  anisotropic 
flow.  It  is  clear  that  the  LaPlace  equation,  (39a),  for  the  field 
between  two  parallel  planes  merely  implies  that 

F = t E (45) 

m 

where  E is  the  applied  electric  field,  independent  of  the  coordinate  x. 
With  Eq.  (45)  and  the  two  functional  relationships,  Eqs.  (41)  and  (43) 
just  derived,  one  is  able  to  replace  n,  0i  and  F of  Eq.  (37)  in  terms  of 
the  mean  velocity  v and  the  applied  electric  field  E.  Equation  (37) 
then  reads 

v 4i  _ A_  4 i _ aE  = 0 (46) 

dx  v3  dx  m 

where 

3k  0i  (o)  [v(o)] 2 

A = 

m 

Since  ©i(o)  and  v(o)  are  the  boundary  values  of  the  temperature  com- 
ponent and  the  mean  velocity  of  the  flow,  A is  independent  of  the 
coordinate  x.  After  the  equation  is  integrated  with  respect  to  x,  we 
have 

V2  + 4 - = G (47) 

v m 
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In  this  equation,  G is  an  integration  constant  independent  of  x.  Since 


at  x = o 


one  obtains  from  Eq.  (47)  that 


G - v‘(o)  + 


With  the  constants  A and  G,  determined  as  above,  Eq.  (47)  can  be  written 


v(x) 

v (o ) 


-M°i  A TyIoTT  l = 2eE> 

'(o)]2  \ |_v(x)J  | m [v(o] 


From  the  above  equation,  one  may  calculate  the  mean  velocity  v as  a 
function  of  the  coordinate  x for  an  electric  field  E. 

In  order  to  expose  the  physical  contents  of  Eq . (48),  it  is  first 
necessary  to  show  how  the  boundary  values  in  the  equation  are  de- 
termined. To  do  this,  the  surface  condition  of  both  cathode  and  anode 
has  to  be  found.  Although  our  understanding  of  the  thermionic  emis- 
sion10 has  been  advanced  considerably  in  the  last  decade,  we  are  still 
unable  to  calculate  the  quantities  n,  v and  0i  accurately  at  the  emis- 
sion surface.  Here  we  use  a reasonable  model  from  which  these  values 
can  be  evaluated.  The  accepted  model  is  as  follows.  All  the  electrons 
which  arrive  at  the  anode  are  collected  by  the  anode  without  causing  any 
stray  disturbances.  The  emission  surface  can  be  simulated  by  a uniform 
surface  of  work  function  W with  a velocity  dependent  transmission 
coefficient  D of  the  simple  form,  Fig.  2,  pg.  4.  Based  on  Nottingham’ s 
experimental  work11,  this  is  designed  to  approximate  the  actual  coef- 
ficient by  a simpler  form  in  which  the  normal  energy  e is  of  the  order 
of  half  an  electron  volt.  Then  the  VDF  at  the  emission  surface  becomes 


3 

f(o)  = d exp  [-  eXp  [ - 

u 3 kT 


m(c i 2 + c22  + c 3 2 ) 


•i  > 0 


=0  , ci  < 0 

where  T is  the  temperature  of  the  emission  surface.  With  this  VDF  at 
the  emission  surface,  the  boundary  values  n(o),  v(o)  and  0i(o)  can  be 


. 


Having  evaluated  the  boundary  values  as  above,  one  can  see  that 
Eq.  (48)  is  normalized  by  the  mean  velocity  of  the  electrons  at  the 
emission  surface.  As  compared  with  Eq.  (2)  of  the  single- valued  veloc- 
ity theory,  Eq.  (48)  contains  an  additional  term,  the  second  term  on  the 
left-hand  side  of  Eq.  (48).  Since  the  factor 

3k  0t(o) 

mv2  (o) 

in  this  term  is  usually  of  the  order  of  unity,  this  additional  term  is 
actually  a correction  term  to  calculate  the  mean  velocity  of  the  elec- 
trons close  to  the  emission  surface.  When  the  ratio 

xLal 

v (o) 

becomes  large,  the  velocities  calculated  from  Eq.  (2)  and  Eq.  (48)  arc 
about  the  same.  After  the  mean  velocity  of  the  electrons  in  the  diode 
is  obtained,  the  number  density  and  the  electron  temperature  can  easily 
be  determined  from  the  functional  relationships  derived  before.  In 
short,  from  Eqs . (41)  and  (43), 

n(x)  = n(o)  [’44] 
v(x) 

and 

0i(*)  = e(o)  2 

v(x> 

As  a conclusion,  a typical  example  for  the  anisotropic  flow  may  Le 
cited  here.  Suppose  the  work  function12  of  the  emission  surface  is 
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1.2  volt;  the  values  of  tp  and  X (Fig.  3)  are  calculated  from  Eqs . 
(33)  and  (34)  against  the  cathode  temperature.  In  a case  of  operation 
where 

cathode  temperature  = 700°K 

cathode  anode  spacing  = 0 . 5 cm 

anode  voltage  with  respect  to  the  cathode  = 300  volts 
one  obtains  from  Eq.  (33) 

tp  = 18  * 10-7  seconds, 

but  from  Eq.  (48)  and  Eq.  (35) 

tj  = 1.9  x 10  8 second, 

The  result  that 

tT  < < tp 

agrees  with  the  specification  given  to  this  type  of  flow. 

3.2  Isotropic  Flow 

As  pointed  out  earlier,  the  number  density  n is  generally  high  in 
the  isotropic  flow.  This  may  be  due  to  a weak  applied  field  which  does 
not  draw  all  the  emitted  electrons  away  quickly  enough.  Also,  it  may 
occur  if  the  emission  current  is  too  high  for  the  applied  voltage  on  the 
anode.  In  either  circumstance,  electrons  may  be  accumulated  to  a very 
high  density  in  front  of  the  cathode.  So  the  problem  in  this  type  of 
flow  would  be  to  solve  simultaneously  Eqs.  (36),  (37),  (38),  (39b)  and 
(40). 

The  processes  of  solving  these  equations  are  very  similar  to  those 
used  in  the  case  of  the  anisotropic  flow.  In  the  first  step,  the 
Poisson’s  Eq.  (39b)  can  be  written  as 

dE  = 4rt  is!  (55) 

ax  mv 

after  relation  (41)  between  the  number  density  n and  the  mean  velocity  v 
is  used.  In  the  second  step,  one  may  replace  n and  0^  of  Eq.  (37)  in 
terms  of  the  mean  velocity  v through  the  two  functional  relationships, 
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In  cm 
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Eqs  (41)  and  (44),  derived  before.  Hien, 

V dl  - _B__  dv  _ p = n 

dx  v 5/3  dx  ( 56.) 

where 

B = -k-e-1-(-0)  [v(o)] 2/3 
3 m 

One  may  eliminate  F from  Eqs.  (55)  and  (56)  in  order  to  obtain  an 
equation  for  the  mean  velocity  v for  this  type  of  flow.  Differentiating 
Eq.  (56)  with  respect  to  x and  then  combining  it  wicn  Eq . (55)  gives 
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(dv\2  _ 4rcje2  = 0 

dx  mv 


(57) 


This  equation  does  not  involve  the  independent  variable  x explicitly13 
and  it  may  therefore  be  reduced  to  a first-order  linear  differential 
equation  by  suitable  transformations.  The  solution  (see  Appendix  IV) 
can  be  written  as 


x + L 


+ 

5v8/3 


(1  - — )vdv 
8 / 2 
v ' 


(58) 


This  solution  contains  two  integration  constants,  L and  M,  which  have  to 
be  determined  by  two  boundary  conditions.  It  is  clear  that  the  non- 
constant force  field,  Eq . (39b),  adds  more  difficulties  in  solving  the 
problem  for  this  type  of  flow. 

Again,  we  shall  use  the  reasonable  model  with  which  estimation  of 
the  boundary  value"  can  be  made  for  the  case  of  the  isotropic  flow.  In 
doing  this,  it  is  obvious  that  the  present  method  is  by  no  means  the 
most  accurate,  but  certainly  serves  as  an  illustration  of  how  boundary 
values  in  this  problem  may  be  determined. 

In  the  isotropic  flow  the  expression  for  the  cathode  emission 


j . ISm  (kT)»  klA  - exp  [- 
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remains  the  same;  but  a large  portion  of  the  mitted  current  is 


(59) 

com- 
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pensated  by  the  returning  component  of  the  current  at  the  emission  sur- 
face. Thus  the  actual  current  density,  j,  of  the  isotropic  flow  is  one 
such  that 

J_  = rs  < < 1 (60) 

•1 1 

This  parameter  Rg  is  the  ratio  of  the  number  of  electrons  actually  drawn 
to  the  anode  and  the  number  of  emitted  electrons  per  unit  area  in  unit 
time.  Hence  Rg  can  be  regarded  as  one  of  the  operating  conditions  of 
the  diode.  In  this  situation,  the  electron  gas  at  the  emission  surface 
wouid  not  follow  the  half-maxwellian  VDF,  Eq.  (49),  but  rather  approach 
a nearly  symmetric  maxweliian  VDF.  Diis  nearly  maxweliian  electron  gas 
may  stretch  itself  out  to  a great  extent  away  from  the  cathode.  In  other 
words,  one  can  then  assume  the  electron  gas  at  the  emission  surface  to 
be  in  equilibrium  with  the  emitter, 

©Jo)  = T (61) 

Since  the  electron  gas  at  the  emission  surface  is  nearly  maxweliian,  the 
returning  component  of  current  is 

n<o>[£tP 

L J 

With  the  above  description  of  this  model,  the  boundary  values  of  n,  v 
and  ©i  can  be  determined  at  the  cathode  surface.  They  are: 


n (o)  = jt(  1 - R 


5 kT 


■(c)  = 

1 - Rg  2nm 


0,  (o)  = T ' (64) 

in  which  jt.  T and  Rg  are  considered  as  known  quantities. 

To  evaluate  the  two  constants,  L and  M,  in  Eq . (58),  another 

boundary  value  of  the  mean  velocity  v must  be  determined  for  this 
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problem.  For  this  purpose  one  may  consider  that  the  number  density  of 
the  electron  flow  decreases  as  the  mean  velocity  increases  in  the  flow. 
From  the  expression  for  the  relaxation  time,  it  is  noticed  that  the 
relaxation  time  is  about  inversely  proportional  to  the  number  density  of 
the  electrons  under  consideration.  Because  the  condition  of  this  type 
of  flow  is  specified  through  the  relaxation  time,  it  is  apparent  that 
the  isotropic  flow  may  become  anisotropic  if  the  electrons  in  the  flow 
are  accelerated  to  reach  sufficiently  high  velocities.  In  the  region 
where  the  mean  velocity  reaches  a sufficiently  high  value,  the  percent- 
age difference  of  the  mean  velocities  calculated  from  the  Eqs.  (47)  and 
(58)  for  the  two  types  of  flow  is  much  smaller  than  the  difference  in 
the  region  of  low  mean  velocity.  Because  of  this  fact,  a good  approxi- 
mate value  of  v(x2)  can  be  calculated  from  Eq.  (47)  for  the  anisotropic 
flow,  provided  that  a suitable  distance  x2  is  selected  for  the  calcula- 
tion. In  the  space-charge- 1 imi ted  flow,  this  distance  x2  may  be  longer 
than  the  actual  cathode-anode  spacing  and  is  used  only  for  the  purpose 
of  detei lining  the  second  boundary  value  of  the  mean  velocity  in  the 
flow. 

The  boundary  values  of  the  mean  velocit’,  at  a distance  t2  and  at 
the  cathode  surface  are  necessary  and  sufficient  tc  evaluate  the  con- 
stants, L and  M,  in  Eq,  (58).  After  evaluating  these  constants,  one  can 
calculate  the  mean  velocity  v(x)  throughout  the  isotropic  flow.  Using 
the  same  manner  as  fqr  the  anisotropic  flow,  one  can  also  find  that 


and 


n (x ) = n(o) 


0!  (x)  = 01  (o) 


V (o) 

|^v  (v) 


V (o)  | 9/3 

V (x) 


(65) 

(66) 


As  a conclusion  for  this  section,  some  calculated  results  may  be 
cited  here  Using  the  same  example  as  given  for  the  anisotropic  flow, 
one  may  raise  the  cathode  temperature  to  1100°K.  When  the  anode  voltage 
is  about  300  volts,  the  parameter  Rg  has  the  value  of  0.005.  From  Eqs. 
(58)  and  (35),  the.  mean  transit  time  of  the  flow  is  obtained 


tT 


2 . 4 x 10  second 


for  an  interval  from  v(o)  to  10  v(o).  However,  from  expression  (33), 


-30- 


APPENDIX  I 

DERIVATION  OF  THE  FLOW  EQUATIONS 


In  this  appendix,  the  derivation  of  the  flow  equations,  (26)-(28) 
is  given  in  detail.  This  is  provided  as  a review  for  the  convenience 
of  the  reader. 

Multiplying  each  term  in  the  Boltzmann  Eq.  (6)  by  p and  integrating 
with  respecr,  to  c in  the  velocity  space,  one  obtains 

' — • ^ £ 

f ip  c ’ Vrf  dc  + / p F * Vcf  dc  = / p dc  (22) 

Using  the  rule  of  integration  by  parts,  one  can  write 

I P Cx  gi-  dc  = -9-  I P Clf  dc  - I fCl  dc  = - nCl^ 

ox1  dxx  ox i oxi  oxi 


(I) 


and 

f P z~-  dc  = If  [^f]  dc?dc3  - If  1—^-  dc  = - n 

oc  i , t = _ a,  dc  i dc  j 

(II) 

where  c and  r are  separate  independent  variables  of  f,  and  c is  inde- 
pendent of  r Since  the  usual  hypothesis  is  that  V-'f  tends  to  zero  a? 
Ci  tends  to  infinity  in  either  direction,  the  integrand  in  the  second 
integral  of  Eq . (II)  vanishes  in  general.  If  p denotes  one  of  the 
summational  invariants,  the  right  hand  side  of  Eq . (22)  vanishes  in 
general.  With  expressions  (I)  and  (II)  available,  Eq . (22)  can  be 
written  in  the  following  foim: 

V ° n P c - n -l  c Vi p + F * V i//  / = 0 (III) 

l J 

It  is  more  convenient  to  specify  the  VDF  and  p as  functions  not  of 
c but  u,  where 

, U - C - c 

Hence  the  next  step  is  to  express  Eq.  (Ill)  in  terms  of  the  peculiar 

*See  for  exau*,-:!e  S.  Chapman  and  T.G  Cowling,  Mathematical  Theory  of  Non-Uniform  Oases, 
Reference  6. 
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velocity  u.  In  changing  the  variable,  we  must  replace 


90  _ 9 0 3JL 

9x  i 9u  9x  j 


c • Vri  jj  - 7uy  • ( c • Vr ) 


One  may  note  that  the  mean  values  are,  of  course,  unaffected  by  the 
change  of  variable.  After  the  substitutions  of  these  terms,  Eq.  (Ill) 
becomes  r ^ 

Vr  • n 0(c  + u ) - n <(c  + u)  4 Vr0  - Vu0  • [(c  + u)  • Vr]  c + F • = 0 


Making  use  of  the  identity 


' [u  * Vr]  c = u vu0  : Vrc 


one  can  reduce  Eq . (IV)  to 


r " ( n>/9  + n 0 V • c + Vr  • n 0 u 


The  notation  of  two  dots  used  in  Eq . (V)  is  the  so-called  double  product 
of  two  tensors.  The  significance  of  the  bar  placed  over  any  quantity 
means  the  mean  value  of  that  quantity. 

In  the  following,  we  shall  consider  Eq  (VI)  for  the  three  known 
values  of  summational  invariants.  For  the  case, 

0!  = 1 


if? 


r x 


1 


Finally  Eq.  (VI)  becomes 


For  the  case 


= 0 
V*  = o 


Vr  (n  c)  = 0 

. “4 

02  = m (c  + u) 


(26) 


then 


and 


02  = me 

n02  Vr  • H = mnr.c  • (^?rc) 

n 02 u = kn0 
Vr  02  = m Vr  H 
u Vr  02  = 0 
Vu02  = m 


Vu  02  u - mu  =0 

Hence,  Eq.  (VI)  with  the  aid  of  Eq.  (26),  in  this  case,  can  be  written 


* (kn  6)  - mn  [F  - c ‘ (Vrc)]  = 0 


For  the  case, 
then 


03  - ^ m (u  + c) 


03_2m^2  + ? k0j 


n 0S u ' kn  0 ■ c 


(27) 
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APPENDIX  II 

FUNCTIONAL  RELATIONSHIP  BETWEEN 
THE  ELECTRON  TEMPERATURE  AND  THE  MEAN  VELOCITY 


In  this  appendix,  we  shall  be  concerned  with  the  acquired  knowledge 
of  the  variation  of  the  electron  temperature  in  a rectilinear  flew.  In 
the  text,  rectangular  coordinates  are  chosen  and  the  rectilinear  flow  is 
along  the  first  coordinate,  say  xt.  In  the  case  of  the  "anisotropic 
flow",  we  have 

T—  02  = j4-  ©a  = 0 (40a) 

dx,  d x i 

and  2 . 

0i  v = constant  (43) 


Since  the  rectilinear  flow  is  along  the  coordinate  Xi,  Eq . (40a)  ex- 
presses the  fact  that  the  temperature  components  perpendicular  to  the 
flow  ©2  and  ©a , remain  constant  throughout  the  flow,  whereas  the  temper- 
ature component  along  the  flow  ©i , decreases  with  the  inverse  square  of 
the  mean  velocity  v.  (See  Eq . (43)).  In  the  case  of  the  "isotropic 
flow", 

©i  - ©2  = 03  (40b) 

a nd 

©i3  v9  = constant  (44) 

Since  in  this  case  the  electron  temperature  components  are  equal  to  each 
other  Eqs . (40b)  and  (44)  imply  that  all  three  temperature  components 
will  decrease  with  the  inverse  two- third  power  of  the  mean  velocity. 

In  the  above  equations  the  variation  of  the  electron  temperature 
is  expressed  in  terms  of  the  mean  velocity  in  a rectilinear  flow.  In 
order  to  exhibit  the  functional  relationships  between  the  electron 
temperature  and  the  mean  velocity  it  is  worthwhile  to  plot  the  varia- 
tions of  the  temperature  components  against  the  mean  velocity  on  log-log 
paper.  In  such  a plot  they  are  straight  lines  as  shown  in  Fig.  4. 

In  considering  the  velocity  distribution  of  the  electrons  in  the 
entire  flow  it  is  necessary  to  study  the  variation  of  the  temperature 
components  in  the  velocity  distribution  function 

02  0 3 


<0  = n(-m-);,/2  ( 0 1 62  0 3)  12  exp  [- 

Zk  k 2 k w - 


Electron  Temperature  Components 


. 


r 


If  the  rectilinear  flow  has  a cylindrical  symmetry  temperature  com- 
ponents perpendicular  to  the  flow  are  equal  to  each  other.  In  this  case, 
it  is  convenient  to  speak  in  terms  of  0Z  instead  of  ©1(  and  0r  instead 
of  02  and  03.  After  this  transformation,  the  velocity  distribution 
function  may  be  written 


( w \ 1 

t m \ % 

2 

mcr 

m(c7  - v)2 

bnkerj  ( 

2rtk0zJ  CXP 

2k@r 



2k0z 

One  may  note  that,  if  the  relationship  between  the  mean  velocity 
and  the  coordinate  x along  with  the  flow  is  established,  the  temperature 
components  can  be  expressed  as  a function  of  x in  the  electron  flow  under 
consideration.  If,  for  example,  the  Chi Id-Langmuir  law  is  accepted, 


then  the  electron  temperature  becomes  a defined  function  of  the  distance 
in  the  flow  Replacing  v,  the  mean  velocity,  in  Eas . (43)  and  (44)  by 
the  Child-Langmuir  relation,  one  obtains: 


anisotropic  flow:  ©!  ® x 4i'S 

and 


02  = 03  - constant 


isotropic  flow: 


0!  « x-‘/e  i = 1,  2,  3. 
These  relations  are  plctted  in  Fig.  5 
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APPENDIX  III 

EVALUATION  OF  THE  CONSTANTS  AT  THE  EMISSION  SURFACE 

In  this  appendix,  the  calculations  will  be  carried  out  under  the 
assumption  that  at  the  emission  surface  the  VDF  is 

f ( o ) = D exp  [-  exp  [-  ^-3  , Cj  > 0 


= 0 


!,  < 0 


and  the  transmission  coefficient  D is 


D = 1 


m ci 

“2^T 


Vi  me  i 2 -i.  e £ 
Vi  me  i 2 < e e 


Let  us  first  integrate  the  VDF  with  respect  to  c2l  c3,  i.e., 

, jZ  ',.C  “c'dc*  ■ 2hV  (tf1)  D exp  t-  ep  exp  [-  ”1^1  (VII) 

The  integration  with  respect  to  Ci  is  as  follows: 

( 2 2 2 

r m lnc*  r mCl  i , ^ r“  r mci  i • i trrT t \ 

J -S~r  exp  [-  — — ] dci  + J exp  [-  -=7=;!  dCl  (VIII) 

o i.eo  2kl  „ , 'A  2kT 


-LT  Jc  expf!l!!£i_] 

2e6  ) 1 ^ ' 2kT J 


<^>* 


tZae.)  2 

• m r mCl  u v 

exp["  TkT  dc*' 


(%i> 


r mc  1 -j  , 

exp  [-  w]  uc, 


With  the  transform 


2kT 


40- 


t 


and 


dCX  = (2JiX)  2 Jy 


expression  (VIII)  becomes 


exp  [-  m + LkT)-,-/2  /(kT>  exp  [-  y ] d>  + (2)il) 

(2me5)/2  kT  (2m)/2e6  ° 


/ exp  [-  y2]  dy 


£ y* 


Since  £ is  about  0.5  volt  and  kT  is  about  0.17  electron  volt  even  at 


2000°K,  the  value  of  (j^j)  is  often  much  greater  than  unity.  The  third 


term  has  the  same  value  but  different  sign  as  the  first  term,  so  only 
the  second  term  is  left  to  be  evaluated.  Expression  VIII  can  therefore 
be  approximated  within  an  error  of  three  percent  by 


_kX  (Hkl) 

2ee  2m 


a 


According  to  Eq . (7)  in  the*  text,  one  obtains  from  expressions  (VII)  and 
(VIII) 


3/9 


n(o)  = (2ntoKT)-_ . (k.T)  exp  [_  eW] 

2 h 3 ee  kr 


Second,  J can  be  evaluated  by  first  considering 
'A 

3 2 2 

lie  - m r'  „ CD  , 

c iexp  [-  dc, 


3 2 

r m • iiiCi  I-  mci  n r 

J -A exp  [-  — — ] dci  t J 

0 2ee  2kT 


o _ e % 
(2ae.) 


2kT 


kTc  1 me/. 

„p  t-  ] 


K 


( kT) 2 r mc 1 i 

_ em.  exp  h 

- exp  [-  i 


m 


2 kT 
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-v  -ti 


* ' . vt  **. 


1 - exp  r-  ] 
meo  k i 


According  to  Eq.  (8)  in  the  text,  one  then  gets 

j , 4gn  [ml  u - exp  {-  «£)]  „p  (-  s® 


and 

v ( o ) = 2 (2-kl)  [1  - exp  (-  p-^)] 

Ttm  k 1 

Third,  the  following  integrals  are  considered  in  evaluating  the 
temperature  component  @i  at  the  cathode  surface: 

. •>  - c .X 
{‘^T)  2 

r ^c i / v 2 / 

J -cT-r  lei  - v ) exp  { 

o zeo 

(IX) 


mci 

~2kf 


)d 


00 


Cl 


+ / 
<1 


(c i - v)?  exp  (- 


me  ■ 


2kT 


) d 


Cl 


Again,  with  the  transform 


me  i 2 

= y 

2kT 


and 


A „ 

^ 1 


(mA  dy 
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expression  (IX)  Lecoii.es 


/(,‘t)  y 9 fy  - (-’a-)  1 v]  exp  [-  v2]dy  + / (-k^) 

J eom  2 k I , m 

O (£.6) 

kT 


3/2 


[y  - ^2kT^  v]2  exP  t”  ys3  dy 


Since  the  value  of  ) is  often  much  greater  than  unity,  the  above  ex- 

k r 


pression  may  be  approximated  by 


(kl)  £kl)3/  ‘ / y9  [y  - (-®-)72  vl  exp  [-  y dy 

o Z k 1 


'A 


(X) 


e u m 


Using  the  known  general  integrals 


j”  yY  exp  [-  a y2]  dy  = - ^ 1 * | * | 


•Lrl8 


-(Y  + 1 )/2 


Y = even 

integer 


r°°  Y r 2,  , I “(Y  + D/2  Y _ 1.1 

J y exp  L-  a y ] dy  = £ a (Y-  „ 1 ) ; 

o Z Z 


Y = odd  integer 


and  the  expression  of  the  mean  velocity  v evaluated  above,  Eq.  (X),  is 
reduced  to 

U 


r3n  -_£)  (kl,  (2iI)J/a 


Qv'n 


e e m 


According  to  the  definition  of  0,  Eq.  (11),  and  the  number  density  r. 
evaluated  above,  one  can  simplify  the  expression  so  that 


0!  - 2.IL_-r  J.  T 

7t 
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APPENDIX  IV 

SOLUTION  OF  THE  EQUATION 


(1 


E— W d— £ + 

v®/3'  dx2 


(1  + -5J3-)  Of*)9  - 

3U8/3  dx  mv 


0 


The  equation  does  not  involve  the  independent  variable  x explicitly 
and  it  may  therefore  be  reduced  to  a first-order  differential  equation 
by  letting 


d x.  = 

dx 


y 


Then 


and  = y 4* 

dx2 


(1  - — E^vy  iX  + (1  + -5_B_)  y3  _ 4TjeI 

8/S  fly  0 8 /3 

v uv  3u  inv 


0 


Furthermore,  let 

2 

y = * 

then  the  above  equation  may  be  written  as 


di 
d v 


2 

+ - 


v 


1 

V 


-JU 

n / 3 


/ 


2 i ) 

mv  j 1 - 
V 


This  is  a first-order  linear  differential  equation,  so  it  can  be  solved 
by  the  general  method.  First,  evaluate  the  integral 


It  is  more  convenient,  if  one  uses 


,9/3  = 


w a 


nd  d\ 


= A 
8 


-s  /a 


dm 


(XI) 
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After  this  substitution,  integral  (XI)  becomes 

J + 5.  B , 
r 3 3d  a 


4 u - B a 


= 3.  f da  + OR  f dJ^ 

4 u a(a  - B) 


= 3-  In  u + 2 In  -B 
4 a 


. i„ 

10/3 


Second,  evaluate  the  integral 


/ _l5ii! (v8/3  - to 9 dv 

„ n 10/3 

mv(1~-^T)  V 

0 / n 


n . 2 t 8 / 3 r»\8/3 

jlius.  / Ly - jl.v — dv 

m „WS 


8^1  [/  dv  - / B~isq 

m ..  8 I 3 


(XII) 
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The  solution  of  the  linear  equation  can  be  written  as 
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1 O f 3 o • 2 OD 

M * + 8JUe_lv  n + 3B  .j 

(v8/3-B)9  m 5va/3 


10/3 

[ V 

(v8/3  - B) 


-1 


10/3 


(v3/3  - B) : 


[M  + ^e— v-  (1 


m 


5v 


aB-)] 

B/  3 


where  M is  an  integration  constant.  So 


dx.  = ± — [M  + 8.5ie-V(i  + _3-B.-)]/2 


dx 


(v 


8/  3 


B) 


c 8/3  ' 
5v  ' 


If  the  positive  sign  is  chosen,  then 


dx  = y.5/3-  [M  + 8-’lje2v (1  + 


d v y »/  3 _ 


5v 


9/3 


and 


-A 


x + L = / [M  * (L  + -££_))  2 {l  - _EL_)  V du 

m ,-8/3  9/3 

m 5v  v • 


where  L is  another  integration  constant.  It  is  clear  that  two  boundary 
conditions  are  necessary  in  order  to  evaluate  the  integration  c ons  tants  , 
M and  L 
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DETERMINATION  OK  ELECTRON  TEMPERATURE 


INTRODUCTION 
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1 . 

Ever  since  the  discovery  of  thermionic  emission,  one-dimensional 
electron  flow  has  been  an  interesting  problem.  The  theory  presented  in 
Part  I of  this  report  deals  with  the  electron  gas  in  a rectilinear 
fiow.  In  the  theory,  the  electron  flow  between  the  cathode  and  the 
anode  can  be  described  completely  with  four  macroscopic  quantities: 

1.  External  force  field  F 

2.  Electron  number  density  n 

3.  Mean  velocity  c 

4.  Electron  temperature  0 

The  first  three  quantities  are  familiar  to  us  since  they  have  been  con- 
stantly used  in  the  s ingle - va 1 ued  velocity  theories.  Experimental 
methods  to  determine  these  macroscopic  quantities  were  developed  long 
ago  and  are  not  discussed  in  this  paper.  However,  the. application  of 
the  last  quantity,  in  this  context  the  electron  temperature  tensor  0, 
may  be  a constructive  contribution  in  this  field. 

In  Part  II  of  this  report,  we  shall  devote  ourselves  to  the  devel- 
opment of  a method  which  would  enable  us  to  measure  the  electron  temper- 
ature in  a rectilinear  flow.  Much  of  the  attention  will  be  given  to  the 
design  of  an  apparatus  with  which  the  experiment  is  carried  out.  The 
present  work  is  Ly  no  means  a complete  study  of  the  subject,  but  the 
results  of  this  experiment  proved  to  be  an  interesting  exposition 
concerning  the  temperature  of  the  electron  gas  in  a parallel  plane  diode 


f 
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2 THEORETICAL  CONSIDERATIONS 

In  order  to  achie.e  the  objectives  stated  in  the  introduction, 
theoretical  consideration  must  first  be  given  tc  the  design  of  an 
apparatus  which  would  be  appropriate  for  the  experiment  to  be  performed. 
The  apparatus  should  have  an  electron  source,  say  a cathode  and  an 
electrode  acting  as  an  anode  to  collect  the  arriving  electrons  Since 
the  electron  flow  discussed  in  Part  I of  this  report  is  of  rectilinear 
type,  the  structure  of  the  apparatus  is  certainly  expected  to  provide 
a paral lei -plane  flow  between  the  cathode  and  the  anode.  In  studying  the 
velocity  distribution  of  the  electrons  in  the  paral le 1 -plane  flow,  use 
can  be  made  of  a field-free  region.  Ideally,  one  would  take  a sample  of 
the  electrons  from  the  para  1 le  1 ••  plane  flow  and  set  this  sample  of  elec- 
trons free  in  a field-free  region.  The  space  distribution  of  the 
current  in  the  field-free  region  would  then  exhibit  the  velocity  dis- 
tribution that  the  electrons  had  before  they  were  set  free.  With  the 
given  definition  of  the  electron  temperature,  an  observed  velocity 
spread  can  readily  be  expressed  in  terms  of  an  electron  temperature. 
This  scheme  would  avoid  the  use  of  a probe  usually  applied  in  obtaining 
information  at  the  point  of  interest.  Here,  in  contrast  to  the  probe 
measurements,  the  point  of  interest  is  placed  at  the  boundary,  hence 
there  is  no  probe  disturbance  in  the  measurements.  The  most  appropriate 
choice  of  a boundary  would  be  the  collector  electrode.  Thus,  allowing  a 
small  sample  of  the  flow  to  pass  through  an  aperture  into  a field-free 
region  would  accomplish  the  desired  function  With  an  arrangement  of 
this  sort  in  mind  the  essential  feature  of  an  apparatus  can  be  sketched 
(Fig  1)  The  apparatus  consists  of  two  parts,  a para  1 le 1 - plane  diode 
and  a field-free  anode  chamber.  Although  the  detailed  description  of 
the  apparatus  will  be  given  in  the  next  section,  its  essential  struc 
ture  can  be  easily  recognized  A small  section  of  all  the  electrons 
forced  in  a Pierce  arrangement  to  flow  parallel  from  the  cathode  to  the 
anode  is  permitted  to  enter  a field- free  region  through  a fine  hole  in 
the  anode  Due  to  their  radial  velocity  component  the  thin  beam  will 
expand  and  the  angular  current-density  distribution  can  be  measured  by 
again  sampling  a small  section  of  that  beam  through  a fine  hole  in  a 
slide.  As  a result  the  current -density  distrilution  at  the  back  wall 
of  the  anode  chamber  gives  us  the  information  about  the  velocity  dis- 
tribution of  these  electrons 


In  this  experimental  arrangement,  deviations  from  the  ideal  case 
described  above  can  be.  expected.  They  will  be  taken  up  in  detail  in 
a later  section  and  will  turn  out  to  be  very  small  indeed.  In  this 
section,  however,  we  shall  restrict  ourselves  to  the  calculation  of  the 
current- density  distribution  at  the  back  wall  of  the  anode  chamber  for 
the  ideal  case. 

For  this  calculation,  Fig.  2 is  an  illustrative  diagram  of  the  beam 
expansion  due  only  to  electron  temperature.  Anode  plane  is  at  z = 0; 
anode  orifice  is  set  at  the  origin.  The  back  wall  of  the  anode  chamber 
is  at  2.  = Z 
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FIGURE  2 DEAM  EXPANSION  DUE  TO  ELECTRON  TEMPERATURE 

The  calculation  of  the  current- density  distribution  at  the  back 
wall  of  the  anode  chamber  is  based  on  the  velocity  distribution  function 
dealt  with  in  Part  I of  this  report  From  the  point  of  view  of  the 
particle  dynamics,  the  momentum  and  energy  of  the  particles  are  the 
primary  variables  in  the  dynamic  problems.  Consequently,  in  practice 
one  would  not  meet  functions  of  the  particle  velocity  higher  than  the 
quadratic  term.  As  shown  in  Part  I of  this  report,  the  moments  higher 
than  the  third  moment  of  the  velocity  distribution  function  do  not 
appear  in  all  the  equations.  Owing  to  this  fact,  one  is  able  to  select 
a representative  VDF  such  as  the  one  given  in  Eq.  (30),  Part  I,  of  this 
report  It  is  noted  that  the  given  function  is  only  a representative 
one  ir.  the  class  of  the.  VDF’s  and  will  approach  the  "true  function'1  up 
to  the  third  moment.  However,  the  VDF  given  in  Eq . (30)  is  the  appro- 
priate one  for  the  purpose  and  can  be  widely  used  in  this  field.  Since 
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the  structure  of  the  apparatus  will  have  a cylindrical  symmetry,  the  VDt 
in  Eq  (30)  is  now  expressed  with  variables  in  the  cylindrical  co- 
ordinates, namely, 


f = n 


/ m \ 

2rtk0r' 


(— m_ 
2"k©z 


% 

) 


exp 


r mcr 

2k6r 


m (cz  - v) ' 

-ife; — 1 


(1) 


where  n is  the  number  density,  v is  the  mean  velocity,  0r  and  ©z  are 
the  temperature  components  in  the  r and  z directions  of  the  electrons  at 
the  anode  orifice. 

With  the  velocity  distribution  function  of  Eq . (1),  the  calculation 
can  be  made  in  a straight-forward  manner.  Per  unit  area  of  the  anode 
orifice,  the  current  between  the  concentric  cones  with  apex  angles  20 
and  2(0  +6)  as  shown  in  Fig.  2 is 


j(M>  ■ /“  dc,  / 
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2rtk0r  2nk0z  0 


c 2 t a r,  6 
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• exp 


mCj.2  m(cz  - v)‘ 

2k0r  2k  0Z 


The  above  integration  over  Cz  is  carried  from  zero  to  infinity  since  only 
the  electrons  with  a position  velocity  component  in  the  z-direction  will 
contribute  to  the  current.  Since  the  area  of  the  anode  orifice  is  Ttr02, 
the  current  density  of  the  Z- plane  between  the  polar  angles  0 and  0+6 
becomes 


.1(6,6) 


j(9,6)Tiro2 

nZ2[tan*(G  + 6 )—  tana0] 


Since  the  angle  6 is  small,  Cr  may  then  be  approximated  by  Cz  tan  0 
under  the  radical  sign.  Ordinarily, 
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go  the  above  expression  can  be  written  as 

To 
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where  I0  is  the  total  current  through  the  anode  orifice, 


Io 


2 

nv  n r0 


Because 


1 - exp' 


lim 

6-0 


tan26)  1 - (1  + ^ tan 2 ( 0 + 6))  ']) 
2k0z  0r  0r 


tan2(0  + 6)  - tan 2 9 


mv 

2k0 


the  expression  of  the  current  density  at  the  Z-plane  with  a polar  angle 
0 can  he  reduced  to 


J (0) 


In  0 „ -3/2  2 I 2 0 „ 

2. (1  + 2 tan  6}  exp  <-  - (1  + ^ tan‘0) 
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For  small  values  of  9,  this  expression  can  be  simplified  to 


J(0)  = — (-n«e_l)  exp 
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This  current -density  distribution  at  the  Z-plane  against  the  polar  angle 
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is  generally  known  as  Normal  or  Gaussian  law.  It  can  also  be  written 
as : 

iog.l-Uil]  ■ - „e!  (4) 

where 

J(0)  = ~ T\ 
rr.Z 2 

and 

2 

m v 

T\  = 

2k0r 

In  this  equation,  logeJ(0)  is  a linear  function  of  02 , and  the  slope  of 
the  straight  line  is  T\ . Since  £q . (2)  is  valid  in  practice,  it  can  he 

found  from  the  flow  equations  in  Part  I of  this  report  that 

% mv?  = eV  (5) 

Therefore,  the  r-component  electron  temperature  can  be  expressed  as 


(6) 


It  may  be  noted  that  the  beam  radius  at  the  Z-plane  is  large  com- 
pared with  the  radius  of  the  anode  orifice.  For  high  velocity  beams, 
this  Z-plane  must  be  at  a great  distance  from  the  anode.  In  this 
experimental  set-up, 

Z = 1 . 200  inch 
2ro  = 0. 0135  inch 

then  Z/2r  = 89.  In  this  case  the  diameter  at  the  Z-plane  is 

o 

expected  to  be  much  larger  than  the  diameter  of  the  anode  orifice.  This 
justifies  the  consideration  made  in  the  calculation  that  the  anode 
orifice  can  be  regarded  as  a point  source  in  this  experimental  arrange- 
ment. If  the  beam  diameter  at  the  Z-piane  is  not  much  larger  than  the 
diameter  of  the  anode  orifice,  a correction  has  to  be  made  for  the 
finite  size-  cf  the  anode  orifice. 


3 . EXPERIMENTAL  APPARATUS 


Ihe  diagram  and  the  apparatus  are  shown  in  3a  .and  >b.  Eusentially, 
the  apparatus,  is  a parallel  plane  diode  (Fig.  1)  which  has  cylindrical 
symmetry  with  respect  to  the  z-axis.  In  this  design,  the  anode  is 
stationary  and  the  cathode  is  movable.  Anode,  detector  and  a toroidal 
ring  are  all  at  the  same  potential  V with  respect  to  the  cathode.  A 
small  orifice  of  0.0135  in.  diameter  is  provided  at  the  center  of  the 
anode.  The  chamber  behind  the  anole  is  almost  closed  and  can  be  con- 
sidered as  a field-free  region.  Embedded  in  the  back  wall  of  the 
chamber  is  a running  slide  on  which  there  is  another  small  orifice  of 
0.0135  inch  diameter.  It  is  this  orifice  which  acts  as  a movable 
opening  on  the  back  wall  of  the  chamber. 

A great  effort  has  been  made  to  have  a parallel  electron  flow  in 
the  diode.  When  the  ca thode -anode  spacing  is  small  compared  to  the 
radius  of  the  cathode,  the  electron  flow  close  to  the  z-axis  can  be 
considered  as  parallel  to  the  z-axis.  It  may  be  noted  that  the  beam 
forming  electrode  is  of  the  Pierce  type,  which  is  helpful  in  maintaining 
a parallel  electron  flow  when  the  cathode -anode  spacing  and  the  poten- 
tial difference  V are  increased.  The  potential  on  the  beam  forming 
electrode  is  so  adjusted  that  the  current  to  the  toroidal  ring  is  small 
or  zero.  The  electron  emitter  is  made  heavier  than  the  ordinary  cathode 
sleeves  to  insure  its  flatness  and  to  increase  the  heat  capacity.  An 
emission  paste  of  half  BaC03  and  half  Sr003  is  used  and  the  top  of  the 
emitter  is  electrophoretically  coated  with  4.4  mg/cm2.  The  coated  sur- 
face looks  very  homogeneous  under  the  microscope. 

Electrons  which  have  emerged  from  the  first  orifice  and  also  get 
through  the  second  orifice  will  then  be  collected  by  the  detector. 
The  disturbances  caused  by  the  anode  orifice  will  be  considered  in 
Section  4.  as  aberrations.  When  the  disturbance  is  small,  the  electrons 
emerging  from  the  anode  orifice  actually  constitute  a sample  of  the 
electrons  arriving  at  the  anode  plane.  In  Section  4,  it  will  also  be 
shown  that  the  flight  of  the  electrons  in  the  chamber  is  affected  very 
little  by  the  mutual  repulsion  force  of  the  electrons.  Therefore,  the 
electrons  having  emerged  from  the  anode  orifice  can  be  considered  as 
flying  through  the  chamber  with  their  initial  velocities  at  the  anode. 

As  a result,  the  current  density  distribution  on  the  back  wall  of 
the  c imber  is  a direct  consequence  of  the  velocity  distribution  of 
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electrons  at  the  ancde.  The  detector  current  indicates  the  current 
density  at  that  point  on  the  back  wail  of  the  chamber,  at  which  the 
orifice  of  the  running  slide  is  set.  By  pulling  the  running  slide  one 
may  obtain  the  current  density  on  the  back  wall  as  a function  of  r,  the 
radial  distance  from  the  z-axis.  From  this  information  the  electron 
temperature  at  the  anode  plane  is  evaluated. 

Since  the  anode  is  supposed  to  collect  all  electrons  which  arrive 
at  the  anode  without  disturbing  the  electron  flow,  the  electron  flow  in 
the  various  cathode-anode  spacings  can  be  considered  t he  same  as  long  as 
the  same  anode  current  is  maintained  by  adjusting  the  anode  voltage. 
When  the  same  anode  current  is  maintained,  the  anode  plane  can  be  placed 
at  any  distance  from  the  cathode  surface  to  determine  the  electron 
temperature  at  that  plane.  For  one  value  of  anode  current,  the  electron 
temperature  can  then  be  determined  as  a function  of  the  distance  from 
the  cathode  surface. 


ESTIMATION  OF  ABERRATIONS 
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As  described  in  the  preceding  section,  the  apparatus  is  designed 
for  the  purpose  of  approaching  the  ideal  case  of  t.  e experiment.  How 
good  an  approximation  is  attained  can  be  seen  in  the  study  of  the 
aberrations  in  this  arrangement.  It  is  conceivable  that  three  kinds  of 
aberrations  would  exist  in  the  experiment.  FirsL,  the  presence  of  the. 
anode  orifice  will  cause  an  electron  icr.s  effect  on  the  emerging  Leatn; 
second,  the  emerging  beam  will  undergo  an  expansion  due  to  its  space 
charge  effect;  third,  the  secondary  electrons  from  the  anode  orifice 
will  also  give  rise  to  a current-density  distribution  at  the  back  wall 
of  the  anode  chamber.  Each  of  these  effects  will  be  estimated  separately 
in  the  following  three  paragraphs. 

4 I Electron  Lens  Effect 

The  estimation  of  the  divergent  action  of  the  anode  orifice  is  made 
as  if  tiiere  were  no  velocity  spread  in  the  electron  flow.  Consider  a 
cylinder  (Fig.  4),  symmetric  with  respect  to  2-axis,  of  cross-section 
area  rtr02  and  length  li  + 12,  just  long  enough  so  that  the  electric 
fields  across  the  ends  of  the  cylinder  are  substantially  uniform  over 
the  area.  Thus,  Ly  Gauss’  Theorem, 

1 2 0 ^ 

2nv0  I Er(zi  dz  = - Ttro^\Ei  - E2  ) - 4~2rc*'  (it  + l2)ne 

-it 

o r 

I Er(z)  dz  = - — - (Ei  - E2)  - 2nr0(li  + l2)ne  (7) 

-It  2 

where  n is  the  number  density  of  electrons  in  the  cylinder,  E;  and  E2 
are  the  electric  fields  at  the  left  and  right  ends  of  the  cylinder,  and 
Er  is  the  radial  component  of  the  electric  field  across  the  mantel  sur- 
face of  the  cylinder. 

Since  the  axial  velocity  of  the  electrons  can  be  regarded  as 
constant  in  the  cylinder,  the  edge  electrons  will  gain  a radial  com- 
ponent of  velocity  equal  to 

“•See  preceding  page. 
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In  view  of  Eq.  (7),  the  gain  of  the  radial  component  of  velocity  is  due 
to  the  additive  result  of  two  divergent  actions  one  comes  from  the 
change  of  electric  field  caused  by  the  anode  orifice,  the  other  comes 
x r cm  tiic  presence  of  space  c ha rgc.  The  first  action  will  be  treated  as 
an  electron  lens  effect14;  the  second  action  as  a space  charge  effect. 

To  estimate  the  electron  lens  effect,  o.ie  must  obtain  an  expression 
for  the  radial  velocity.  From  Eq.  (8) 


ei'o(Ei  - E2 ) 


2mv 


(9) 


More  explicitly,  each  electron  entering  the  left  end  of  the  cylinder  has 
the  same  velocity  as  given  in  Eq . (5)  in  the  direction  of  the  z-axis. 
Passing  through  the  region  of  the  cylinder,  each  electron  will  gain  a 


radial  component  velocity  as  given  in  Eq . (9).  Then  their  ratio  will 


bi 


~v~ 


r0(Ei  - E2 ) 


4V 


where  Ei,  in  a space  charge  limited  diode,  is 


E 


4Y 

3D 


and  E2 , in  the  anode  chamber,  is 

E2 


0 


To  describe  the  electron  lens  effect,  we  may  best  express  the 
effect  in  terms  of  a polar  angle  0i  (at  the  center  of  the  orifice)  which 
the  cross-sectional  area  of  the  beam  at  the  back  wall  would  subtend  In 
doing  so,  sc  have 
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In  Fig.  6,  this  polar  angle  9a  is  plotted  against  —-.  It  should  be 

r0 

noted  that  the  electron  path  depends  only  upon  the  shape  of  the  poten- 
tial field  and  not  upon  the  magnitude  of  the  potential.  As  long  as  the 
configuration  is  maintained,  the  change  of  the  anode  potential  should 
not  cause  any  change  of  the  deflection  angle. 

^*2  Space  Charge  Effect 

Based  on  the  analysis  given  by  Spangenberg1 5,  one  can  also  estimate 
the  beam  spread  due  only  to  space  charge  effect.  For  the  purpose  of 
analysis,  it  is  convenient  to  start  at  a point  where  the  radial  velocity 
of  the  electrons  can  be  considered  zero.  Assumed  conditions  are: 

a.  Electrons  in  the  beam  are  distributed  in  a rotational  symmetry 
with  respect  to  the  z-axis, 

L.  The  axial  velocity  of  the  electrons  in  the  beam  remains  the 
same  in  the  consideration. 

The  analysis  was  primarily  developed  for  high  velocity  electron  beams. 
If  there  were  no  velocity  spread  in  the  electron  flow,  the  beam  spread 
thus  predicted  for  the  case  of  low  velocity  electron  beams  should  still 
be  reasonably  accurate. 

When  the  diode  (Fig.  5)  is  space-charge- limited , the  following 
relation  exists: 

» / 
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Z , and 


is  its  radius  at  the  orifice  Although  Eq.  (11)  cannot  be  expressed  in 
terms  of  elementary  functions,  it  is  clear  that  the  percentage  increase 
in  the  beam  radius  depends  only  upon  the  ratio  of  Z and  D.  Since  r0  and 
Z are  known  in  this  experiment-  . Eq . (11),  for  the  convenience  of  indica- 
ting the  functional  relationship,  can  also  be.  written 

Y = Y ( 0/ro ) 
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Putting  the  space  charge  effect  on  the  same  basis  for  comparison, 
we  have 

r riY 

p«,  - arc  tan  — 

it 

In  Fig.  6,  this  polar  angle  05  is  also  plotted  against  D/ro-  If  the 
diode  is  in  space- charge- limited  operation,  the  beam  expansion  due  to 
the  space  charge  effect  i?  also  independent  of  the  anode  potential. 

From  Eqs.  (7),  (8)  and  (11),  the  radial  velocity  that  an  electron 
would  gain  before  reaching  the  back  wall  is  the  sum  of  the  velocities 
due  to  the  electron  lens  effect  and  those  due  to  the  sp  ce  charge  effect. 
In  Fig.  6,  the  sum  of  the  polar  angles,  0i  and  02-,  is  also  plotted 
against  D/r0.  Because 

tan  6i  ' tar.  02  - tan  (9i  + 02),  (0i  + 02)  < 90° 


(Oi  + 02)  is  a maximum  polar  angle  estimated  for  the  beam  expansion  ir 
this  experiment. 

1.3  Secondary  Electron  Effect 

In  order  to  examine  the  effect  of  secondary  electrons  produced  at 
the  anode  orifice,  we  shall  calculate  the  current- densi ty  distribution 
at  the  back  wall  of  the  anode  chamber  due  to  the  secondary  electrons 
alone.  In  this  calculation,  the  emitted  secondaries  are  considered  to 
enter  the  field-free  region  with  negligible  mutual  interactions. 
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Figure  7 diagram  of  the  anode  orifice 

The  average  nnmher  nf  pl.ctrons  per  second  that  will  hit  the  inner  nail 
of  the  anode  orifice  (Fig.  7)  is: 
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Taking  the  yield  factor  to  be  unity,  secondary  electrons  thus  produced 
can  be  expressed  as: 
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(12) 


In  the  work  of  secondary  electron  emission,  it  has  been  reasonably 
accepted  that  the  angular  distribution  of  secondary  electrons  follows 
the  cosine  law,  independent  of  the  angle  of  incidence  of  the  primary 
electrons.  In  other  words,  the  number  of  secondaries  emitted  per  unit 
solid  angle  is  greatest  in  the  direction  normal  to  the  emitting  surface 
and  decreases  with  increasing  angle  of  emergence  0 as  cos  9.  Consider- 
ing the  secondary  electron  er  ission  per  unit  area  of  a plane  surface  at 


FIGURE  8 SECONDARY  ELECTRON  EMISSION  FROM  AN  OR  I FI CF 
the  origin  (Fig.  8)  one  can  write 
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in  the  primed  coordinates,  or 


= h cos  9 


(13) 
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in  the  unprimed  coordinates,  where  is  is  the  secondary  electron  current 
per  unit  area  of  emitting  surface,  Q is  the  solid  angle  which  P subtends 
at  the  origin,  and  h is  a proportional  constant. 

To  describe  the  secondary  eLectron  emission  from  the  inner  wall  of 
an  orifice,  we  consider  that  the  current  at  P is  the  sum  of  the  contri- 
butions by  a small  plane  surface  in  one  revolution  with  respect  to  the 
axis  of  the  orifice.  ’jfliis  is  equivalent  to  the  case  where  the  emitting 
surface  stays  stationary  but  P is  revolving  around  the  Z-axis.  Then  the 
current  at  P would  be  the  same  as  the  current  to  the  ring  (generated  by 
P in  this  manner)  from  the  small  plane  surface  at  the  origin.  Thus,  Eq. 
(13)  is  written: 
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This  describes  the  angle  dependence  of  the  secondary  electron  emission 
per  unit  area  of  a cylindrical  surface  at  the  origin. 

Now,  the  secondary  electron  emission  of  an  orifice  can  Le  expressed 

£»  5 


H sin2  0 
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where  Ig  is  the  secondary 
proportional  constant.  By 


electron  emission  of  an  orifice,  and  K is  a 
normalizing  the  expression,  one  finds  that 
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In  deriving  Eq.  (14).  it  is  assumed  that  tw  <<  r0  <<:  B.  Let  Js  be  the 
current  density  at  the  back  wall  of  the  anode  chamber  due  to  the  second- 
ary electrons  alone,  then 
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From  Eq.  (140,  we  have 
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then 
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This  gives  the  current  density  at  the  back  wall  of  the  anode  chamber  as 
a function  of  the  polar  angle  at  the  center  of  the  anode  orifice,  if 
only  secondary  electron  emission  processes  are  considered.  In  Fig.  9, 
the  value  of  sin2  9 cos 3 6 is  plotted  against  0. 

In  order  to  see  the  effect  of  the  secondary  electrons  in  our 
measurements  of  current  density  distribution  at  the  Z-piane,  Eq.  (15) 
can  be  written  as 
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in  the  iight  of  Eq.  (12).  In  this  experiment,  for  example, 
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and 


r\  = 26.6 


then 


J-S-(  9 ).  i 3 x 10 
J(0) 


This  ratio  is  very  small  as  compared  with  that  in  Eq.  (4).  It  will  be 
seen  that  the  ratio,  J(0)/J(O),  encountered  in  the  measurements,  is 
always  greater  than  10-1.  Hence,  the  secondary  electron  effect  is 
entirely  negligible  in  this  experiment. 

4.4  Conclusion 

In  the  study  of  this  arrangement,  estimations  of  the  aberrations 
have  been  made  for  each  of  these  effects.  It  is  rather  fortunate  that 
the  beam  expansion  predicted  from  the  existing  electron  temperature 
alone  is  much  larger  than  all  the  effects  arising  from  the  aberrations. 
This  will  also  be  seen  in  the  next  section  where  the  experimental  re- 
sults are  presented.  Since  the  current-density  distribution  at  the  back 
wall  of  the  anode  chamber  is  predominantly  due  to  the  electron  tempera- 
ture, the  velocity  spread  of  the  emerging  electrons  can  readily  be 
expressed  in  terms  of  an  electron  temperature. 

4.5  Experimental  Results 

Before  we  present  the  current-density  measurements  at  the  back  wall 
of  the  anode  chamber,  it  would  be  proper  to  assure  that  the  assumptions 
made  in  the  theoretical  considerations  are  fulfilled.  Particular  atten- 
tion has  been  paid  to  the  activation  of  the  large  surface  cathode,  and 
the  proper  performance  of  the  diode  has  been  carefully  checked.  Measure 
ments  of  the  space  charge  flow  for  different  cathode-anode  spacings  have 
been  carried  ou.t.  The  results  are  shown  in  Figs.  10  and  11.  The 
straight  lines  in  booh  figures  represent  the  theoretical  Chi Id -Langmuir 
space-charge  law  in  logarithmic  coordinates;  first,  in  the  form 
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FIGURE  10  REDUCED  CURRENT  DENSITY  AS  A FUNCTION  OF  CATHODE  ANODE  SPACING 
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second,  in  the  form 

jAD:  = QV3//  2 

where 

jA  is  the  current  density  in  pa/cm* 

V is  the  anode  potential  in  volts, 

D is  the  cathode-anode  spacing  in  cm, 

Q is  a constant,  2.335. 

The  points  in  both  figures  represent  the  observed  values  reduced  to 
the  corresponding  units.  It  may  be  noted  that  the  micrometer  drive  has 
a zero  reading  of  about  0.018  cm.  This  would  cause  the  last  two  sets  of 
points  in  Fig.  10  to  deviate  from  a straight  line,  but  the  observed 
values  indicate  very  strongly  that  the  Child-Langmuir  space-charge  law 
is  definitely  applicable.  Since  in  Sections  4.1  and  4.2,  the  Child- 
Langmuir  law  has  been  used  in  the  deductions,  it  is  important  to  verify 
that  the  diode  operates  properly  under  appropriate  conditions. 

Having  established  the  proper  functioning  of  the  diode  with  respect 
to  its  main  characteristics , one  may  now  proceed  to  the  determination  of 
the  desired  electron  temperature.  As  was  shown  in  Section  2 of  this 
Part  (Eqs.  (4)  and  (6)),  the  electron  temperature  at  a given  point  of 
interest  is  evaluated  from  the  investigation  of  the  current  distribution 
on  the  back  wall  of  the  anode  chamber.  The  corresponding  equations  were 

1 °Se  [j-J-g-]  = --  ".e2  (4) 

k0_  = & (6) 

r n 

where  .J(fl)/,I(0)  is  the  relative  current  density  nn  the  Lack  wall  of  the 
anode  chamber,  0 is  the  polar  angle,  r\  is  the  slope  of  Eq . (4),  V is 
the  anode  potential  and  ®r  is  the  r-component  electron  temperature. 

In  order  to  obtain  the  desired  data,  the  cathode  is  set  at  a 
definite  temperature;  with  the  cathode  micrometer  drive,  the  cathode  is 
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set  at  a desired  spacing;  and  a stable  anode  current  is  adjusted  by 
arranging  for  an  appropriate  anode  potential.  The  current  distribution 
along  the  back  wall  of  the  anode  chamber  is  then  determined  by  measuring 
the  current  passing  to  the  collector  for  different  positions,  of  the  fine 
aperture  on  the  running  slide.  With  the  given  geometry  of  the  anode 
chamber  it  is  easy  to  express  the  position  of  the  aperture  on  the 
running  slide  in  polar  angles  9 with  respect  to  the  anode  orifice.  To 
facilitate  a comparison  between  the  experimental  data  and  the  theoreti- 
cal expressions,  Eqs.  (4)  and  (6),  the  measured  values  must  be  expressed 
in  terms  of  the  logarithm  of  the  relative  current  density,  J(9)/J(0), 
and  the  square  of  the  polar  angle,  9'. 

A large  set  of  measurements  have  been  carried  out  and  some  of 
the  typical  results  are  shown  in  Fig.  12.  The  points  in  Fig,  12  repre- 
sent. the  measured  values  snd  show  c good  indication  of  constituting  a 
straight  line  as  required  by  Eq.  (4).  These  results  seem  to  indicate 
very  strongly  that  the  spread  of  the  electron  beam  behind  the  anode 
orifice  is  caused  by  the  radial  velocity  distribution  of  the  electrons 
at  the  moment  of  entering  the  anode  chamber,  rather  than  by  the  dif- 
ferent aberrations  as  lens  action  of  the  anode  orifice,  space  charge 
effect,  and  secondary  emission.  For,  if  there  were  no  velocity  spread 
in  the  flow.,  the  beam  cross  - sec  tion  at  the  back  wall  of  the  anode 
chamber  would  (a)  subtend  a polar  angle  less  than  (9t  + 02  ) at  the 
center  of  the  anode  orifice,  and  (b)  be  independent  of  the  anode  poten- 
tial applied.  But  nothing  of  that  sort  coulu  be  observed.  First,  the 
beam  cross-section  at  the  back  wall  of  the  anode  chamber  subtends  by  a 
far  greater  polar  angle  than  (9-  + 02)  at  the  center  of  the  anode 
orifice.  Second,  the  beam  expansion  has  a strong  dependence  on  the 
anode  potential.  This  makes  the  velocity  spread  of  the  electrons  the 
most  important  factor  in  accounting  for  the  beam  expansion,  end  confirms 
the  estimated  small  values  of  those  aberrations  as  given  in  Section  4. 

It  was  mentioned  in  Section  3 that  in  spite  of  different  spacings 
between  cathode  and  anode,  the  same  regions  of  the  electron  flow  may  be 
considered  as  having  the  same  characteristics,  provided  the  cathode 
temperature  and  the  current  density  in  the  diode  are  the  same.  There- 
fore, this  method  allows  measuremen l of  the  electron  temperature  at  any 
arbitrary  point  within  a parallel  plane  flow  which  is  defined  by  just 
two  parameters,  namely,  cathode  temperature  and  current  density.  In 
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order  to  obtain  the  desired  electron  temperature  at  a given  point  of 
interest  along  the  £iow  which  is  characterized  by  its  cathode  tempera- 
ture and  its  current  density,  one  has  to  place  the  anode  plane  at  the 
point  of  interest  and  determine  the  electron  temperature  by  observing 
the  current  distribution  or.  the  Lack  wall  of  the  anode  chamber  as 
described  above.  Varying  the  cathode-anode  spacing  and  adjusting  the 
anode  potential  such  that  the  anode  current,  remains  constant  gives  a 
complete  description  of  the  electron  temperature  as  a function  of  the 
distance  from  the  cathode.  It  may  be  noted  that  the  so  determined 
electron  temperature  at  each  point  of  interest  is  an  evaluation  of  a set 
of  data,  obtained  by  measuring  the  current  distribution  at  the  back  wall 
of  the  anode  chamber  for  the  point  in  question,  as  described  above. 

A set  of  experiments  was  carried  out  to  establish  the  relationship 
between  the  electron  temperature  and  the  distance  from  the  cathode  sur- 
face for  a given  cathode  temperature  and  anode  current.  The  results  of 
these  experiments  for  a constant  cathode  temperature  of  1275°K,  and  for 
different  current  densities  in  the  flow  are  represented  in  Tig.  13.  One 
interesting  fact  may  be  immediately  noted,  namely,  the  decrease  of  the 
electron  temperature  with  the  increase  of  the  distance  from  the  cathode 
surface.  This  decrease  of  the  temperature  can  easily  he  accounted  for 
if  one  remembers  that  in  the  theoretical  part  of  this  report,  a decrease 
of  the  temperature  with  increasing  distance  from  the  cathode  surface  was 
predicted  for  the  isotropic  flow  where  the  three  temperature  components 
arc  equal  to  each  other.  The  corresponding  equation  reads: 


fir  ® 


(1.44) 


Although  this  equation  establishes  only  a correlation  between  the  elec 
tron  temperature  and  the  mean  velocity  of  the  electron  stream  at  the 
point  of  consideration,  it  is  easy  to  convert  the  expression,  with 
t he  aid  of  the  Child- Langmuir  equation,  into  a temperature- distance 
relation  In  Appendix  II  of  Part  I this  relationship  is  derived  and  ac- 
cording to  Eq.  2,  Part  II,  one  obtains  for  the  isotropic  flow  the  desired 
correlation  between  the  electron  temperature  and  the  distance  from  the 
cathode  surface, 
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Since  the  data  in  Fig.  13  are  represented  in  log-log  coordinates, 
any  exponential  dependency  would  present  itself  as  a straight  line,  the 
slope  of  which  would  correspond  to  the  exponent  of  the  independent 
variable  (see  also  Part  1,  App.  II).  For  comparison,  a -4/9  slope 
representing  the  isotropic  case  is  indicated  in  the  same  graph.  Al- 
though for  large  distances  there  is  an  indication  that  the  electron 
temperature  may  follow  a -4/9  power  law,  for  smaller  distances  this  is 
obviously  not  the  case.  At  this  point  no  account  can  be  given  for  the 
deviation  from  the  expected  dependency.  Only  a much  more  extensive 
experimental  investigation  of  the  region  close  to  the  cathode  surface 
can  provide  sufficient  data  to  allow  further  theoretical  considerations. 
They  will  be  presented  iri  a later  report. 

It  may  be  interesting  to  compare  the  experimental  results  obtained 
in  this  study  with  the  well  known  fact  that  the  electron,  temperature 
should  increase  with  increasing  plat.e  voltage.  Since  in  the  usual 
experimental  arrangement  the  cat hnde -anode  spacing  is  kept  constant,  the 
above  statement  is  equivalent  to  the  statement  that  the  electron  temper- 
ature should  increase  with  increasing  current  density.  Figure  14  is 
obtained  by  arranging  the  experimental  data  such  that  the  electron 
temperature  is  plotted  against  the  current  density  with  the  distance  as 
parameter.  One  may  note  that  the  electron  temperature  tends  to  reach  an 
asymptotic  value  for  high  current  densities. 

One  of  the  most  startling  findings  in  the  experimental  results  of 
this  study  is  the  extraordinarily  high  electron  temperature  as  compared 
to  the  cathode  temperature.  To  make  this  point  particularly  clear,  in 
Fig.  15  the  electron  temperature  at  a constant  distance  of  0.13  ’em  is 
plotted  against  the  cathode  temperature  for  two  different  plate  voltages, 
50  and  100  volts.  Although  the  cathode  temperature  does  not  exceed 
1300°K,  the  temperature  of  the  electron  gas  at  the  anode  orifice  assumes 
values  of  the  order  of  10,000  to  40,000°K.  Before  giving  an  explanation 
of  this  phenomenon  we  should  note  that  an  agreement  of  the  temperature 
of  the  electron  gas  with  the  temperature  of  the  cathode  is  usually  found 
by  the  retarding  field  method  at  low  cathode  temperatures.  In  the 
evaluation  of  the  data  obtained  from  the  retarding  field  measurements, 
an  important  assumption  is  made,  namely,  thcit  the  space  charge  in  the 
diode  is  negligible  Due  to  this  assumption,  the  range  of  these  meas- 
urements is  usually  limited  to  low  cathode  temperatures.  In  fact  at 
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FIGURE  lb  RELATION  BETWEEN  ELECTRON  TEMPERATURE  AND  CATHODE  TEMPERATURE 


high  cathode  temperatures  the  experimental  curves  do  not  show  tiie 
expected  abrupt  bend,  hence  the  evaluation  of  the  electron  temperature 
from  these  curves  seems  rather  doubtful.  Since  the  electron  temperature 
from  high  temperature  cathode  has  not  beer,  experimentally  determined, 
the  above  results  do  not  contradict  the  agreement  iovre  at  tne  iow 
cathode  temperatures. 

Before  concluding  the  experimental  part,  an  attempt  will  be  made  to 
explain  the  extremely  high  temperature  of  the  electron  gas  close  to 
the  cathode.  This  explanation  will  rest  essentially  on  the  experimental 
evidence  that  the  temperature  of  the  electron  gas  decreases  sharply  if 
the  temperature  of  the  cathode  is  lowered,  as  indicated  in  Fig.  15. 
Although  the  present  apparatus  is  not  sensitive  enough  to  measure  the 
electro/;  temperature  for  very  low  cathode  temperatures,  it  seems  reason- 
able to  believe  that  the  temperature  of  the  electron  ge.s  will  approxi- 
mate the  cathode  temperature  in  the  lower  regions  This  assumption  may 
be  well  supported  by  a closer  examination  of  the  emission  mechanism  of 
oxide-coated  cathodes. 

It  is  known  that  the  oxide  coating  prepared  in  the  conventional 
processes  has  a porous  structure  the  porosity  ranging  from  65%  to  about 
85%12.  Experimental  evidence,*6  that  the  outer  layers  of  the  oxide 
coating  are  chiefly  responsible  for  the  electron  emission,  was  found. 
Figure  16  may  represent  a schematic  diagram  of  such  an  oxide - coa ted 
surface.  This  structure  suggests  that  at  low  temperature  the  surface 
grains  may  possibly  be  responsible  for  the  entire  emission.  The  re- 
placement of  the  electrons  lost  by  the  surface  grains  through  the 
process  of  emission  will  be  supplied  by  electrons  coming  from  the  deeper 
layers  through  a process  of  conduction  in  the  oxide  coating.  Since  the 
electron  emission  in  this  case  may  be  considered  as  originating  from  the 
grains  at  the  surface,  the  velocity  spread  or  the  temperature  of  the 
emitted  electrons,  would  correspond  to  the  thermal  temperature  of  the 
oxide  coating  as  in  the  case  of  electron  emission  from  pure  metals 

Experiments ^ have  shown  that  conduction  through  the  oxide  coating 
is  due  to  two  mechanisms  acting  in  parallel- -the  electronic  conduction 
through  the  grains,  predominating  at  low  temperatur es  and  the  conduc- 
tion through  the  electron  gas  in  the  pores  between  the  grains  which  is 
preponderant  at  high  temperatures  Since  at  higher  temperatures  the 
conductivity  of  the  coating  will  decrease,  on  the  other  hand  the  emis- 
sion will  rapidly  increase  with  increased  temperature  :i.t  is  clear  that 
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at  elevated  temperatures  the  conduction  will  lag  behind  the  emission. 
With  this  picture  in  mind  it  is  not  difficult  to  see  that  at  high 
temperatures  a part  of  the  emission  may  come  from  the  electron  gas  in 
the  pores  between  the  grains.  In  the  process  of  emerging  from  the 
narrow  channels  between  the  grains,  the  electrons  will  probably  suffer  a 
series  of  collisions  with  the  heavy  grain  particles  and  may  assume  tem- 
peratures which  by  far  exceed  the  temperature  of  the  grain  particles.  18 
It  is  interesting  to  note  that  the  above  explanation  can  also 
account,  for  the  fact  that  the  electron  temperature  will  increase  if  a 
higher  anode  potential  is  applied.  In  a recent  paper19  by  Loosjes,  Vink 
and  Jansen,  it  was  found  that  under  pulsed  operation  the  velocity  spread 
of  the  emitted  electrons  from  oxide  coating  would  often  amount  to 
hundreds  of  electron  volts. 


Oxide  Coating 


Vacuum 


FIGURE 


16  ILLUSTRATIVE  DIAGRAM  OF  THE  OXIDE  COATING 
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